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PREFACE TO PART HI. 

The third part of the work, now offered to the pro- 
fession, contains problems in thermodynamics of per- 
fect gases and saturated steam. The first part con- 
tains problems in machine design, the second part 
problems in hydrauUcs and the remaining two parts 
are devoted to mechanics of materials and electrical 
engineering respectively. Thus in using this work a 
student or an engineer who wishes to review calculus 
or analytics, or to acquire facility in applications of 
higher mathematics to engineering problems, may 
select at first the part of the work which deals with 
problems in that branch of engineering with which he 
is most familiar, or in. which he is particularly interested. 

The book as a whole is not intended to bring out 
anything new, either in mathematical methods or in 
practical deductions. The author's aim was simply 
to collect and to arrange in a systematic way the va- 
rious applications of analytic geometry and of calculus, 
already in use. The book may be called a summary 
of the most conamon engineering appUcations of higher 
mathematics, or a mathematical cross-index to engineer- 
ing text-books. It fulfills its purpose if it saves the 
teacher the trouble of consulting many engineering 
books for the purpose of selecting a few mathematical 
problems for his students. The author also hopes that 
the book may stimulate interest in higher mathematics 
among his fellow engineers and thus help to a better 
understanding of some intricate relations where at 
present ^^ rule of thumb " prevails in design. 

• • • 
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iv PREFACE. 

The reader will find the author's views on teaching 
mathematics to engineering students in the preface to 
Part I and in the Dialogue following that preface. 
He is also referred to Part I for a list of reference works 
on mathematics and for an Appendix entitled '^What 
a Senior in Engineering ought to know about Mathe- 
matics. " 

The author wishes to acknowledge gratefully the 
assistance of his former students, Mr. J. G. Pertsch, 
M.E., Instructor in Electrical Engineering at Cornell 
University, who read the manuscript on thermody- 
namics, and Mr. Fred G. Switzer, M.M.E., who solved 
most of the problems. Mr. A. C. Stevens, M.E., In- 
structor in Cornell University, read the proofs, and made 
a number of valuable suggestions. 

Cornell University, Ithaca, N. Y. 
January^ 1916. 
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ENGINEERING MATHEMATICS. 



CHAPTER I. 
WORK DONE BY EXPANDING A PERFECT GAS. 

(a) Fundamental Cases. 

The action of the two most important classes of 
prime movers, viz., steam engines and gas engines, is 
based upon the expansion of a gas or a vapor. In the 
steam engme, the working gas is water and vapor, in 
the gas engine it is mostly atmospheric air brought to 
a high temperature and pressure by an admixtm^e of 
certain explosive gases or vapors. In order therefore 
to understand the action of these prime movers and 
be able to calculate their dunensions for the perform- 
ance of certain work, it is necessary to know some of 
the properties of gases. 

To make the conditions more definite it is customary 

to consider a imit weight of gas, say a kilogram (or a 

pound). Then the properties of an arbitrary quantity 

of the same gas may be deduced by simply multiplying 

or dividing the results obtained for the unit weight by 

the number of units of weight contained in the gas under 

consideration. Thus, for instance, if the total mass 

of gas contains G units of weight, then the total volimie 

V = Gv\ total heat H = Gh; total work W =^ Gw; etc., 

where v, h, w, etc., denote the respective quantities 

3 



4 ENGINEERING MATHEMATICS. [Chap. I. 

per unit we^t of gas. For our purposes, a kilogram 
(or pound) of gas is characterized by three quantities: 

(1) volume; 

(2) pressure; 

(3) temperature. 

This is illustrated in Fig. 1: Let one kilogram of a 
certain gas, for instance, atmospheric air, be comprised 
in a cylinder C provided with a piston B, pressed down 
by a total force of P kilograms.* 
If the cross section of the cylin- 
der, or the area of the piston, 
f is S square meters, then P/S = 

\ p represents the pressure of the 

i gas in kilograms per square 
meter. Let the cylinder be 
surrounded by a bath so that 
the gas within is maintained at 
a constant temperature t. Then 
the gas will assume a definite 
a volume V. In this sense, it is 
said that a gas is perfectly de- 
fined by two out of its three characteristics, say, by its 
pressure and temperature. 

Before taking up changes in a gas when two of its 
properties, for instance volume and temperature," are 
varied at will, let us investigate what happens to the 
gas when one of its three characteristics is kept constant, 
while one of the other two properties follows variations 
dependent upon changes made at will in the remaining 
characteristic. It will be understood that we shall 
have in mind only perfect gases, that is, such that are 

• The absolute pressure P includea the force of atmospheric pressure 
and the weight of the piston. 




Chap. I.] 



EXPANDING A PERFECT GAS. 



far from their point of liquefaction and which therefore 
obey the laws of Boyle and of Gay-Lussac. 

Changes at Constant Volume. 

Let the temperature of the gas (Fig. 1) be gradually 
increased by supplying heat from outside; the gas will 
expand, raising the piston, and in order to keep the 
voliune constant it will be necessary to increase the 
force P. Another way to perform the experiment is 
to lock the piston in place and to measure the pressure 
with a pressure-gage, such as is used on steam boilers. 
Plotting results as shown in Fig. 2, a straight line is 

p 
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Temperature 



Fig. 2. — Variation of pressure with temperature, at a constant 

volume. 

obtained, indicating how the pressure increases with the 
temperature. Producing the lines to their intersec- 
tion with the ^-axis, a point is obtained, called the 
absolute zero, or a temperature at which the gas exerts 
no pressure whatever. Temperatures counted from 
this point are called absolute temperatures, and are 
denoted by T, in contradistinction to ordinary Centi- 
grade temperatures denoted by U Experiments on so- 
called perfect or permanent gases give practically the 
same value of the absolute zero temperature, namely 
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273° C. below that of freezing water. Thus, for in- 
stance, if the temperature t = 25° C, then T = 25 + 
273 = 298° C. abs. Fig. 2 shows that, at a constant 
volume, the pressure in a perfect gas is proportional 
to the absolute temperature, or 

P = C^T, (1) 

where Ci depends on the kind of gas and on the imit 
of volume selected. Experiment shows also that the 
amoimt of heat h conmiimicated to the gas per unit of 
its weight is proportional to the rise in temperature, or 

h = c.{U-ti), (2) 

where Cv is the so-called specific heat at constant vol- 
ume. It will be seen from Eq. (2) that Cv, by definition, is 
the amount of heat necessary to raise the temperature 
of one kilogram of gas one degree Centigrade, the vol- 
ume of the gas remaining constant. 

In the metric system the temperatures in Eq. (2) 
are supposed to be on the Centigrade scale, while heat 
content h is measured in calories per kilogram weight 
of gas, one calorie being the amount of heat necessary 
to raise the temperature of one kilogram mass of 
water one degree Centigrade, at the temperature of its 
highest density (between 4 and 5° C). The specific 
heat of atmospheric air at constant volume is 0.169. 
This means that only 16.9 per cent as much. heat is 
necessary to raise the temperature of air as is required 
to raise the temperature of an equal weight of water 
an equal nmnber of degrees. 

The volimae of the vessel being constant, no external 
work is done by the piston, and all the heat appUed is 
used to raise the temperature of the gas. According 
to the modem views on the constitution of gases, par- 
ticles of a gas are in constant motion in all directions 
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(kinetic theory of gases), and an increase in temperature 
means an increase in the mean velocity of the par- 
ticles. It is said therefore that the heat suppUed in 
this case mcreases the intrinsic or internal energy of 
the gas without producing external work. 

Prob, I. — How much heat is required to raise the temperature 
of 3 kg. of air 85° C. at a constant volume? 

Ans. 43.095 kg-cal. 

Changes at Constant Pressure. 

Let now the piston (Fig, 1) be loaded with a con- 
stant weight Pj and let heat be communicated to the 
gas from the outside; the piston will travel upward, 
and at the same time the temperature of the gas will 
increase. The heat suppUed is partly converted into 
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Fig. 3. — Variation of volume with temperature at a constant 

pressure. 

the external work of raising the weight, and partly in- 
creases the temperature or the intrinsic energy of the 
gas. Plotting volimies against temperatures, a straight 
line is obtained (Fig. 3) which intersects the temper- 
atiu'e axis at a point corresponding to the absolute zero. 
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the same as in Fig. 2. Thus, analogously to Eq. (1) 

we have 

t; = C2r, (3) 

where the coeflScient of proportionaUty C2 depends on 
the kind of gas and on the pressure used (Gay-Lussac's 
law). 
The heat communicated may be represented here as 

h =Cp(fe-<i), (4) 

where Cp is the specific heat of the gas at constant 
pressure. It is evident that Cp > c», because in this 
case the heat not only raises the temperature of the gas, 
but also performs some external work. Let the piston 
rise by an amoimt I; the total external work performed 
is then PI = pSl = p (v2 — Vi) kg-m. 

If volume and pressure are used as abscissae and ordi- 
nates (Fig. 4) the change in the state of the gas is repre- 




Volume 



Fig. 4. — Change in volume at a constant pressure. 

sented by a horizontal line, which shows that the pres- 
sure of the gas is independent of its volume. The work 
p {v2 -- Vi) done by the gas is in this case proportional 
to the area of the rectangle viV2ba, since its base is 
(v2 — Vi) and the altitude represents the pressiu'e p. 
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Chap. I.] EXPANDING A PERFECT GAS. 9 

According to a fundamental law of nature, heat is 
equivalent to work; that is, in order to perform a cer- 
tain amount of work, a proportional amount of heat 
must be used. In the metric system this relation is 

1 large calorie = 427 kg-m. 

It is customary to denote the inverse value of this 
ratio by A, so that in the metric system 

A = 1/427. 

Thus, the amount of heat per imit weight of gas, used 
in raising the weight P, is 

Ap {v2 — Vi) heat units. 

On the other hand, the amount of heat necessary to 
raise the temperature of a perfect gas from <i to fe is 

Cv (k - tl)j 

independent of whether the gas is at a constant volume 
or not.* Hence, the total heat given to the gas is 

h = Cp(t2- ti) = Cvik- ti) +Ap(v2- vi), . (5) 

where Cp is known as the specific heat at constant 
pressure. 

Prob. 2. — 3 kg. of air at a pressure of 7 metric atmospheres (1 
metric atmosphere = 1 kg. per sq. cm.) and at a temperature of 
15° C. occupies 0.362 cu. m. The air is allowed to expand at 
the constant pressure until its volume becomes 0.563 cu. m. It is 
found that its temperature is 175° C. From these data find how 
many calories were added to the gas, and what percentage of the 
heat was used in performing the external work and in increasing 

* This is true of perfect gases only; with imperfect gases part of the 
heat serves to disintegrate particles of gas and neither raises its tem- 
perature nor performs external work. This is called the internal latent 
heat of the gas. For perfect gases the latent heat is zero, as has been 
proved by the famous "porous plug" experiment of Joule and Kelvin. 
Hence, no matter what change of state the gas is undergoing, the in- 
crease of its intrinsic energy always corresponds to the increase in 
temperature, and does not depend upon the volume or pressure. 
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the intrinsic energy of the air. Also calculate Cp, knowing that for 
air Ct, = 0.169 (in the meter-kilogram system). 
Am. 114.06 kg-cal.; 28.88 per cent; 71.12 per cent; Cp=0.2375. 
Prob. 3. — In the preceding example, the weight was lifted uni- 
formly in 12 sec. What horse-power was the engine developing 
during the process? One metric horse-power is equal to 75 kg-m. 
per sec. Ans. 54 metric h.p. 

Changes at Constant Temperature. 

Let the cyUnder (Fig. 1) be a perfect conductor of 
heat and be surrounded by a very large bath maintained 
at a constant temperature, so that in whichever way 
the pressure and volmne of the gas may change, its 
temperature remains constant. Let the gas be in equi- 
librium occupying a volume t; at a pressure p. Remove 
a smaU part of the Weight P; the piston will be raised 
upward by a small amount and an equilibrium estab- 
lished again. At the same time a small quantity of heat 
will flow into the cylinder from the outside to supply 
the work of lifting the weight. On the other hand, if 
the weight on the piston is increased, the gas will be 
compressed, and the work of the weight will be con- 
verted into heat. Since the gas is surrounded by a 
bath maintained at a constant temperature, this extra 
heat will flow into the surrounding medium. These 
changes in volume with pressure at a constant temper- 
atm-e are called isothermal changes, and are repre- 
sented by curves in Fig. 5. Each curve refers to 
a particular temperature. Experiment shows that 
under such conditions pressm*es are inversely propor- 
tional to volumes. 

Let the initial and the final pressures and volumes of 
the gas be pi, p2 and Vi, V2 respectively; let some inter- 
mediate values be p and v. Then 

PiVi = pv = P2V2 = const. . ... (6) 
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The value of the constant depends upon the kind of 
gas and the temperature at which the experiment is 
performed (Boyle's law). Eq. (6) represents analyti- 



CONST. 




yolame 
FiQ. 5. — Changes in a gas at a constant temperature. 

cally a family of equilateral hyperbolae, shown in Fig. 5, 
each curve corresponding to a certain temperature. 

To calculate the work performed in this case it is 
necessary to integrate over infinitesimal displacements 
of the piston, because the pressure is diflferent in every 
position. The work done during an infinitesimal dis- 
placement dl is 

P •dl = pS •dl = p •dv, 

so that the total work performed, per imit weight of 
gas, is 

w = I p *dv. 

This work is represented graphically in the p-v diagram 
(Fig. 5) by the area of the figure ViV2ba. This is the 
reason why pressure and volume are usually selected 



w 
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as coordinates in considering changes in state of sub- 
stances, since the external work performed is repre- 
sented by the area of the diagram between the hmits of 
volume. 

As the expansion is isothermal and no heat is neces- 
sary in perfect gases to disintegrate particles, all the 
heat supplied is converted into work, and we have, in 
heat units, 

h = Aw = A j pdv. 

Substituting the value of p from Eq. (6) gives 

= h/A = piVi I dv/v = piVi Ln (f;2M).* . (7) 

From this equation the heat suppUed and the work per- 
formed may be calculated when the initial conditions 
of the gas and the Umits of expansion are given. If, 
instead of V2 and vi, the limits of pressure pi and p2 are 
given, the equation becomes 

w = h/A = pit;iLn (P1/P2). ... (8) 

Prob. 4. — With the initial data given in Prob. 2, assume that 
the air expands isothermally instead of at a constant pressure. 
What will be the work performed and the heat supplied; also what 
is the final pressure of the gas? At what volume and pressure 
should the expansion begin, in order to get the same mechanical 
effect as in the case of expansion at a constant pressure of Prob. 2? 
Ans. 33,600 kg-m.; 78.6 kg-cal.; 4.5 metric atmos.; 

pi = 12.36 metric atmos.; Vi = 0.2047 cu. m. 

Prob. 5. How many kilogram-meters of work are required to 
compress isothermally 3 kg. of air which originally occupies 425 
liters at the temperature of freezing water, to double the original 
pressure? One kilogram of air occupies a volume of 0.7733 cu. m. 
at atmospheric pressure and 0° C. Atmospheric pressure = 1.0333 
metric atmos. = 1.0333 kg. per sq. cm. Ans. 16,500 kg-m. 

* The symbol Ln is used here to denote the Naperian logarithm; 
the common logarithm is denoted by the symbol log. 
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Characteristic Equation of Perfect Gases. 

In the preceding processes one of the three charac- 
teristics of the gas is considered constant, and the 
experimental results are expressed by the three equa- 
tions, (1), (3), and (6). These equations may be 
combined mto 

pv = RT, (9) 

which is called the characteristic equation of a per- 
fect gas. When the volume is constant, the pressure 
increases as the absolute temperature, as required by 
Eq. (1); when the pressure is constant, the volume in- 
creases proportionately to the absolute temperature, as 
required by Eq. (3) ; when the temperature is constant, 
the product pv is also constant, in accordance with Eq. 
(6). 72 is a physical constant which depends on the 
kind of gas. Its value may be determined from Eq. (9) 
if a set of corresponding values of p, v, and T is known. 
For instance, for air, according to Prob. 5, we have at 
the atmospheric pressure and the temperature of freez- 
ing water: Vq = 0.7733 cu. m. per kg.; To = + 273 = 
273° C. abs.; po = 1.0333 X 10^ = 10,333 kg. per sq. 
m. Hence, substituting into (9), 

R = 10,333 X 0.7733/273 = 29.27. 

Eq. (9), being a relation between three variables, can- 
not be represented by a curve, but is represented by 
a surface shown in Fig. 6. Cross sections by planes 
corresponding to t; = const, give a series of straight 
lines which represent changes at constant volume, as 
in Fig. 2. Similarly, cross sections by planes perpen- 
dicular to the p-axis give expansions at constant pres- 
sure, according to Fig. 3. Intersections with planes 
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perpendicular to the T-axis give equilateral hyper- 
bolae of isothermal changes, as in Fig. 5. 

An excellent idea of the properties of a perfect gas 
is obtained by actually making a model of the above 




Fig. 6. — A surface showing the relationship between pressure, volume, 

and temperature in a perfect gas. 

surface, according to Eq. (9). Planes of coordinates 
can be made out of wood; isothermal curves are made 
in the form of ribs out of paste-board, wood, or alumi- 
num; and strings or wires are stretched to represent 
j>-T and v-T lines.* 
A satisfactory picture of the properties of a perfect 

* Such a model has been constructed by the author. See Sibley 
Journal of Engineering, Vol. 29 (1914), p. 12. 
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gas is obtained by projecting isothermals on the p-v 
plane, as in Fig. 6, so that for any p and v, the corre- 
sponding temperature, T, may be read directly, or 
estimated. 

Prob. 6. — Show by reasoning that the values of R for different 
gases are inversely as their densities. 

Prob. 7. — Prove from Eqs. (5) and (9) that 

Cp = Cv + AR. . . . . . . . (10) 

Adiabatic Expansion. 

Besides the processes considered above, when one of 
the quantities p, v, or T remains constant, there are 
innmnerable other ways of changing the condition of 
a gas. Among these processes, one that is of particu- 
lar practical importance is the so-called adiabatic 
change of state, when no heat is communicated to the 
gas from outside. Imagine the cylinder and the piston 
in Fig. 1 to be made of perfect non-conductors of heat ; 
if the pressure on the piston is increased or reduced 
the gas is compressed or expanded. In expanding it 
lifts the weight and performs mechanical work. Since 
no heat is conmiunicated from the outside, this work is 
done at the expense of the intrinsic energy of the gas 
itself; consequently, during an adiabatic expansion the 
gas is cooled. On the contrary, during an adiabatic 
compression some work is performed upon the gas, and 
as the mechanical energy thus converted into heat can- 
not be conducted away, the temperature of the gas 
rises. 

The problem is as follows: being given the original 
pressure, volume, and temperature of the gas, to find 
the law according to which they vary during an adia- 
batic process. As these three quantities are always 
connected by Eq. (9), only two of them need be deter- 
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mined. Since all the three characteristics of the gas 
vary continually, it becomes necessary to consider an 
infinitesimal change. The work p • dv, done during 
such a change, and the consequent decrease of temper- 
ature dT are connected by the relationship 

A'pdv= -c^dT ..... (11) 

since the work is done entirely at the expense of the 
internal energy of the gas. The minus sign before dT is 
necessary because T decreases when v increases. Eq. 
(11) may be obtained directly from Eq. (5) by letting 
fe = and reducing the interval fe — <i to an infinitesimal. 
Eq. (11) is the differential equation of an adiabatic 
change m a perfect gas; the three variables which enter 
into it are connected by the ' characteristic Eq. (9) . 
Hence, eluninating one of the variables between these 
equations, a relation is obtained between the remaining 
two variables. It is easy to eliminate p by dividing 
Eq. (11) by Eq. (9). The result is 

A dv/v = - c,dT/RT, .... (12) 
or, integrating, 

{AR/c) Ln {thM = Ln {T^/T^). 
Taking into account Eq. (10), this is reduced to 

'■^^^)Mw} ■ ■ ■ »« 

The ratio of the specific heats Cp to Cv occurs quite 
often in the theory of gases and is sometunes denoted 
by one letter k, that is, 

Eq. (13) becomes 

(<c - 1) Ln iih/vr) = Ln {T./T^), 
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or 

(v2/vi)^' = Ti/T2 (14) 

This is the required relationship between v and T dur- 
ing an adiabatic change. It may be represented in the 

form 

Tgt^-i = Tivi'^-i = Ttf-"^ = const. . . . (15) 

A relationship between p and v is obtained by eliminat- 
ing T between Eqs. (9) and (15). The result is 

p-iVi" = piVi"^ = pif= const. . . . (16) 

FinaUy, the relationship between p and T is 

y^K/(ic-i) y^-c/cit-i) y«/(K-i) 



P2 Pi P 



= const. • (17) 



The work of an adiabatic expansion is found directly 
by integrating Eq. (11), namely 



p .dv = —Cv j dTj 



or 

Aw = c. (Ti - Ta) (18) 

This result could be easily [foreseen, since the work is 
performed at the expense of the internal energy of the 
gas. Eq. (18) may be expressed through the volume 
and the pressiu'e, instead of the temperature, by using 
Eqs. (9), (15), (16), and (17). 

Prob. 8. — Prove that for an adiabatic process 

Prob. 9. — Knowing that for air, R = 29.27 and c^ = 0.169, 
find Cp and prove that k = 1.405. Ana, Cp = 0.2375. 

Prob. 10. — Demonstrate that on a p-v diagram an adiabatic 
curve (pv* = const.) at any given point is steeper than the isothermal 
(pv = const.) passing through the same point. 
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Piob. II. — Let the gas in Prob. 2 expand adiabatically, instead 
of at a constant pressure. What is the amount of heat converted 
into mechanical work? Ans. 81.12 kg-cal. 

Prob. 12. — On a sheet of crossHsection paper plot p-w curves 
of expansion, (a) at constant pressure, according to Prob. 2; 
(6) isothermal, Prob. 4; (c) adiabatic, Prob. 12. Mark on the 
curves the initial and the final temperatures, the work performed, 
and the heat supplied. 



CHAPTER 11. 

WORK DONE BY EXPANDING A PERFECT GAS, 

(b) Polsrtropic Changes. 

Besides the fundamental cases of expanding a gas, 
considered in the preceding chapter, there are innumer- 
able other combinations of simultaneous changes in 
pressm^e, volume, and temperature. The only condi- 
tion which connects these three quantities is the char- 
acteristic equation of the gas 

Vv^RT, (1) 

so that two out of three of the quantities may be varied 
arbitrarily; the third is then determined from Eq. (1). 
As an example of such thermodynamic processes we 
shall consider in this chapter changes during which 
the quantity of heat supplied is proportional to the 
change in temperature (poly tropic changes). 

During a change of state at a constant volume the 
heat supplied to or withdrawn from the gas is pro- 
portional to the change in temperature, as is expressed 
by Eq. (2) in the preceding chapter. There are an in- 
finite number of other processes in which the volume 
varies during the change and yet the quantity of heat 
supplied is proportional to the change in temperature. 
The problem is to find the general relation between p 
and V for such processes. The condition just stated 
leads to the equation 

dh = cdT = Cvdt + A 'pdv, . . . (2) 

19 
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which indicates that an increase in heat dh is propor- 
tional to the rise in temperature dT, and also that the 
heat suppUed partly increases the intrinsic energy of the 
gas by the amount c« dT, and is partly converted into 
external work p • dv. The coefficient of proportionahty 
c has the dimension of specific heat, and may be called 
the specific heat of a polytropic process, as distinguished 
from the specific heats Cp and Cv at constant pressure 
and at constant volume respectively. 

Since the problem is to determine the relation be- 
tween p and Vy temperature T must be eUminated from 
Eq. (2) by means of the characteristic Eq. (1). By 
dififerentiating the latter we have: 

RdT = pdv + vdp, .... (3) 

which, substituted into Eq. (2), gives 

(Cv — c) {pdv + v dp) + ARp dv = 0. 

But, according to Eq. (10) of the preceding chapter, 

so that substituting and simplifying we find 

^ + n^ = 0, (4) 

p V 

where 

n = ^^^. (5) 

The integral of Eq. (4) is 

Ln (P1/P2) + n Ln (vi/t^z) = 0, 
or 

PiVi"" = P2V2'' = pv" = const. ... (6) 

This is the required equation of a polytropic curve 
between p and v. (Note that the exponent n appUes 
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to V only, and not to the whole expression.) The gen- 
eral shape of such curves is shown in Fig. 7. When 
n = (or c = Cp), p = const., and the curve becomes 
a straight line parallel to the axis of abscissae. For 
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• Constant 
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Vs 
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Vi 

Volume 
Fig. 7. — Polytropic changes in a perfect gas. 

n > the pressure decreases when the volume increases, 
and vice versa. For n < both the volume and the 
pressure decrease and increase simultaneously. When 
n = 1, Eq. (6) becomes identical with Eq. (6) of the 
preceding chapter, and the change of state is isothermal. 
When c = 

n = — = Kj 

and the curve becomes an adiabatic; see Eq. (16) in the 
preceding chapter. 
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Prob. I. — Plot a curve showing the relationship between c and 
n from c=— oo toc = +oo. Explain the physical meaning of 
various values of c, and the corresponding shape of polytropic 
curves. 

Piob. 2. — Prove that for a polytropic curve the ratio of an ab- 
scissa Or (Fig. 8 ) to the corresponding subtangent rg is constant and 
equal to the exponent n of the curve. This is one of the methods 




Fig. 8. — Determination of exponent n. 

for determining this exponent when the curve is given as an experi- 
mental result. 

Prob. 3. — Find the relation between T and v by eliminating p 
from Eq. (2) and integrating the result. Show that the answer is 
the same as when p is eliminated directly between Eqs. (1) and (6). 



Ans. 



Ti ^ (vjy-^ 

T2 \Vi) 



I Prob. 4. — Find the relation between T and p by eliminating v 
from Eq. (2) and integrating the result. Show that the answer is 
the same as when t; is eliminated between Eqs. (1) and (6). 

Construction of Poljrtropic Curves. 

It is sometimes required in the design of a gas engine, 
an air compressor, etc., to construct a polytropic curve 
according to Eq. (6), the exponent n and the initial 
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■ 

pressure and volume being given. This can be done 
by assuming various values of p and calculating the 
corresponding v's with the acid of a table of logarithms, 
or it can be done graphically, using so-called loga- 
rithmic paper (Fig. 9). This is a cross-section paper 
on which divisions are proportional to logarithms of 
numbers, instead of being proportional to the niun- 
bers themselves.* When plotted on such paper aU poly- 
tropic curves become straight lines of different slope. 
This becomes evident by taking logarithms of both 
sides of Eq. (6), viz.: 

log p + n log i; = log pi + n log Viy . . (7) 

or denoting log t; by x and log p by y 

y +nx = const. 

This is the familiar equation of a straight line, in which 
(— n) is the tangent of the angle which the Une makes 
with the X-axis. 

An example will make the method clearer: Let air 
be compressed in the cylinder of an air compressor 
from atmospheric pressure and a temperature of 17° C. 
to a pressure of 6 metric atmospheres. If the com- 
pression were adiabatic the air would be heated con- 
siderably; to reduce the work of compression the air 
is cooled during the compression by spraying water 
into it. Assume that this makes the compression 
follow the law 

pt^'^ = const (8) 

instead of the adiabatic law 

p^.405 = const. 

* If logarithmic croes-section paper is not available, the divisions 
on the upper scale of a slide rule may be used for the scale of coordinates. 
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It ia required to construct the p-v curve according to 
Eq. (8). 

From the characteristic equation we find at the at- 
mospheric pressure (pi = 1.0333 X 10* kg. per sq. m.), 
i^ = 29.27 X (17 + 273)/10,333 = 0.8215 cu. m. 
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v^ volume 111 CU. m. per kg. weight of gas 
Fio. 9. — Use of logarithmic croBa-Bection paper for plotting poly- 

On the curve sheet (Fig. 9) the stretch ab corresponds 
to specific volumes from 0.1 to 1.0 cu. m. per kg. weight 
of gas, and the stretch fee to specific volumes from 1.0 
to 10.0 cu. m. per kg. weight of gas. Along the axis of 
ordinates ad gives pressures from 10* to 10^ kg. per 
sq. m., and de from 10* to 10* kg. per sq. m. The point 
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g is selected to correspond to Vi = 0.8215 and pi = 
10,333. Then the line gh is drawn at a slope of —1.2 
to the axis of abscissae. This line represents the re- 
quired polytropic curve plotted to a logarithmic scale. 
It can either be used directly or re-plotted to an ordi- 
nary scale by simply transferring a sufficient mmaber 
of points. At the pressure of 6 metric atmospheres 
(6 X 10^ kg. per sq. m.), the specific volume, as deter- 
mined from the curve, is 1^2 = 1.90 cu. m. per kg. weight 
of gas (point h). The corresponding temperature is 
found from the characteristic equation: 

7^2 = 1.90 X 6 X 10V29.27 = 389, 

or <2 = 389-273 = 116° C. 

The use of logarithmic paper enables one to judge 
whether or not a given curve of change of state is a 
polytropic curve, and if not, how near it comes to being 
one. This is done by transferring a sufficient number 
of points from the given curve (say, from an indicator 
card of a gas engine or an air compressor) to a 
logarithmic sheet. If all the points he on a straight 
line the curve is of the form pv = const., and the 
inclination of the straight line gives directly the value 
of n. Otherwise, it is possible to tell if n increases or 
decreases with the increase in volume, and also draw 
a straight Une representing a polytropic curve which 
comes sufficiently close to the given curve. 

If a curve is known to be polytropic and it is merely 
required to find n, we have, from Eq. (7), 

^^ Hv^-lo^p. .... (9) 
log t^ - log t;i 

where the subscripts 1 and 2 refer to any two given 
points on the curve. 
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Prob. 5. — In the case shown in Fig. 9 what would be the final 
volume and temperature if the gas were allowed to be compressed 
adiabatically? Same isothermally? Transfer the three curves 
from logarithmic to ordinary cross-section paper in order to become 
familiar with the actual shape of the curves. Why is the work of 
compression greatest for the adiabatic curve and smallest for the 
isothermal curve? 

Am. Adiabatic, U = 209** C; V2 = 0.235 cu.m.; 
Isothermal, <2 = 17** C; V2 = 0.1417 cu. m. 

Prob. 6. — Can the curve gh (Fig. 9) be drawn to a larger scale, 
using the same kind of logarithmic paper? 
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Fig. 10. — A change which is only approximately polytropic. 



Prob. 7. — Determine if the curve ab shown in Fig. 10 is poly- 
tropic, and if not what is the equation of the nearest polytropic 
curve? What error would be made in determining n by directly 
using Eq. (9) between the two extreme points of the curve? 
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Work and Heat Values in Poljrtropic Changes. 

Integrating Eq. (2) we get : 

h=c{T2- Ti) = c. (T2 - Ti) + Aw, 

where h is the total heat per unit weight of gas, suppUed 
or abstracted between the states 1 and 2, and 



w 



= Tpdv (10) 



is the external work done by the gas or upon the gas 
between the same states. We have thus 

h = c(T2-Ti), (11) 

Aw = {c - c«) (7^2 - Ti), . . . (12) 

or both the heat suppUed and the work performed are 
proportional to the change in the temperature of the gas. 
If it is desired to express h and w through pressures 
and volumes, instead of temperatures, Eqs. (1) and (6) 
may be used. Thus, combining Eqs. (1) and (12) gives 

or, remembering that AR = Cp — c», 

^^ c^:^ (^t^__j> , , , (13) 

Cp-Cv (PlVi ) 

In some cases in engineering practice the original 
pressure pi and volume Vi are known, as well as the 
travel of the piston. Therefore the final volume V2 is 
also known, while the final pressure p2 is not inamediately 
given. It is therefore desirable to eliminate it from the 
above equations by means of Eq. (6). Substituting, 
and taking Eq. (5) into account, the result is 



w 



-^U-M"'l . . . (14) 

n — 1 ( \V2/ ) 
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This formula is used to a considerable extent in thermo- 
dynamic calculations. 

When gas is being compressed, instead of expanding, 
Pi and Vi in formula (14) represent the final, and not the 
initial, values (Fig. 8). If it is desired to use the initial 
values p2 and V2 in this case, formula (14) may be 
transformed by means of Eq. (6), which gives pit;i = 
P2V2 (t^/t^i)'*"^ Substituting, we find 



w 



= -2^S(^y"-li. . . . (14a) 
n-ll \vi/ ) 



Prob. 8. — Derive formula (14) directly from Eqs. (10) and (6). 

Prob. 9. — Prove that the ratio of the area piabp2 (Fig. 8) to 
the area Viabv2 is equal to the exponent n of the polytropic curve. 
This gives a method for determining n by means of a planimeter. 
Hint. Area piabp2 = area Viobvz -\- PiVi — p^2' 

Prob.. 10. — When n = 1, expression (14) assumes an indeter- 
minate form 0/0. Evaluate w in this case according to the general 
rule given in differential calculus for indeterminate expressions of 
this form. The result must check with Eq. (7) of the preceding 
chapter, since f orn = 1 the change of state is isothermal. 



Equation (14) in the Form of an Infinite Series. 

Expression (14) is not convenient for numerical 
computations when the value of n is near unity. In 
this case, both the expression in the brackets and the 
denominator tend to become nearly zero, so that a 
small error in the values of Vi and «;2 or of n may 
appreciably change the result. It is therefore 
convenient to expand the expression in the brackets 
into an infinite series of powers of (n — 1). The 
value of (n — 1) being small, the series converges 
rapidly, and a few terms usually suflice in practical 
calculations. 
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The expansion is performed according to the well- 
known Maclaurin series 

fix) =/(0) +f /'(O) +|]/"(0) +|^/"'(0)+ • . . . 
In our case, denoting (n — 1) by x, we have 

/ (0) = 1 ; /' (0) = Ln ?J ; /" (0) = (in JjY; etc., 



V2 \ V2J 

so that 



2 

+ 



\V2/ 1 Vi 2 1 \ Vi' 

In this expression use is made of the transformation 

Ln^=-Ln^; 

Vi Vi 

Vi being, by assmnption, smaller than «^, Ln (vi/th) is a 
negative quantity and may, therefore, for the sake of 
convenience, be replaced by a positive logarithm of 
the inverse ratio. 
Thus, we finally have : 

(15) 

This expression is preferable to Eq. (14) when n is near 
unity, because the terms in the brackets, from the 
second on, are small as compared to 1, and a slight error 
in their calculation is not liable to affect the result 
appreciably. When n = 1 Eq. (15) simply becomes 

w = piViLn—> 

Vi 



• • > • 
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This is another way of evaluating the indeterminate 
form referred to in Prob. 10. 

Prob. II. — Calculate the expression w/piVi, using both formulss 
(14) and (15), in order to judge about their relative convenience and 
accuracy. Assume the ratio of expansion V2 : vi to be equal to 4 : 1, 
and n = 1.05. 

Prob. 12. — It will be seen from Fig. 7 that the work done be- 
tween the volumes vi and V2 is smaller, the larger the exponent n, 
provided that the initial conditions are the same in all cases. Prove 
the same from formula (14) ; in other words, prove that with given 
values of pi, vi, and V2 the value of w increases when n decreases. 

Solution, — It is required to prove that 

an 
The given expression is equivalent to 

ivhere a; = w — 1 and r = -^ > 1. 

Vi 

We have 

dji _ —r' +1 + xliwr 

dx xV 

The denominator is positive; to prove that the numerator is 

always negative let us find its maximimi and minimum values. 

Thus, if 

z = — r* + 1 + xLdt, 

^=Lnr(l-r*) = 0; 

in other words, the critical value of the function is for a? = 0; at 
this value 2 = 0; the second derivative, 

cPz 



dx^ 



= - (Lnr)V, 



is negative, hence the value x = corresponds to a maximum of the 
function z. Thus, the maximum value of 2 is zero, hence z is always 
negative. This proves that 

t«>- 

or that w decreases when n increases. 
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Work Expended in an Air Compressor. 

A single-acting air compressor is represented sche- 
matically in Fig. 11. When the piston P is being moved 
upward, atmospheric air lifts the valve s and fills the 
cylinder. During the downward stroke, valve s is 
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Fig. 11. — Schematic diagram of an air compressor. 



closed and the air is compressed during a part of the 
stroke to the desired pressure pi, after which the valve t 
is opened and the compressed air is discharged into the 
storage tank K in which a constant pressure is main- 
tained by means of the weighted piston Q. Then the 
same operations are repeated. The problem is to 
calculate the mechanical work done by the piston 
during one cycle comprising two strokes. Referring 
to Fig. 8, the net work performed per kilogram of air 
is represented by the area piabjh* During the suction 
stroke, the atmospheric air performs the work OpJ^th] 
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during the compression period the piston perfonns the 
work represented by aftt^vi, and finally during the dis- 
charge period the piston perfonns the work equivalent 
to Opiavi. Thus, the net work is 

w = abv2Vi + Opiavi — Op2&t^ = Pidbp2. . (16) 

The compression curve ab would be an adiabatic if 
the cylinder were impervious to heat; in reaUty the 
curve deviates considerably from an adiabatic, and for 
practical purposes it may be represented by a poly- 
tropic curve in which the exponent n varies for different 
compressors. Thus, using expression (14a) for the 
area abt^Vi, Eq. (16) gives 



w 



= -^S (^V ' - 1 ! +.m - P2V2. . (17) 

n — 1 ( KvJ ) 

In this equation p2 represents the atmospheric pressure, 
Pi the pressure of the compressed air, V2 the volume 
which one kilogram of air occupies at the atmospheric 
pressure p2 and at a certain temperature T2 which must 
be given. Volume Vi corresponds to the pressure pi and 
to a certain temperature Ti which may be calculated 
from Eqs. (1) and (6) if the exponent n is known. 

Piob. 13. — Show that Eq. (17) may be simplified to 

Piob. 14. — Deduce Eq. (18) directly from Eq. (14a) on the basis 
of the theorem stated in Prob. 9. 

Prob. X5. — A double-acting * compressor, working at a speed of 
75 r.p.m., must deliver 1000 kg. of air per hour at a pressure of 

♦ This means that similar valves t and s (Fig. 11) are also placed on 
the top of the cylinder, so that air is compressed during both strokes. 
This doubles the work done with the same siae of cylinder. 



or to 

w 
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10 metric atmospheres (absolute). The air is taken at ordinary 
pressure and at room temperature (say, 20** C). To determine: (a) 
the displacement of the cylinder; (6) kilowatts of the driving 
motor, assuming the efficiency of the compressor to be 75 per cent; 
(c) the temperature at which air is delivered; (d) per cent stroke 
at which the valve t must be opened. The exponent n may be as- 
sumed equal to 1.25. 
Ans. (a) 0.0925 cu. m.; (b) 89.5 kw.; (c) 189° C; (d) 83.73 %. 



Compound Compressors. 

Expression (19) shows that the work necessary for 
compressing one kilogram of air depends not only 
upon the initial and the final pressures, but also upon 
the final volume Vi of the air. Now, with a given 
pressure, the volimae depends on the temperature; 
therefore, in order to reduce the work of compression, 
the temperature during the compression must be kept 
as low as possible. This can be done by dividing the 
process of compression between two cyhnders, and cool- 
ing the air on its way between the cylinders. Thus, 
referring to Fig. 11, the air is compressed in the first 
cylinder to a pressure pm, intermediate between pi 
and jh, and is deUvered to the tank K where it is 
cooled as much as possible at a constant pressure, 
and then is deUvered to a second compression cylinder 
where it is compressed to a higher pressure pi. Such 
an arrangement is called a compoimd or a two-stage 
compressor. 

The problem is to select the intermediate pressure 
Pm so as to reduce the total work of compression to a 
minimum. 

Applying Eq. (19) we get 

(n-l)/n -] r/'n, \(n-l)/n 

It/ \ If iLfm ( 

w 



=n-^S4fe)'-"-]M©' 
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or, introducing the corresponding temperatures T2 and 
Tm which are known, 



w 



=„-^S''{©'""""--']-«"""-']S-<-) 



The variable part of this expression, containing pm, 
is 

/7} \(n-l)/n /-DtX^w-D/w 

and this must be a minimum if the work w is to be a 
minimum. The value of pm which makes y a minimum 
may be f oimd according to the general rule, by equat- 
ing dy/dptn to zero. However, it is much simpler to 
apply the theorem that when the product of two fimc- 
tions of a variable is constant, the sum of these func- 
tions is a minimum when the two fimctions are equal 
to each other (see Prob. 20 below). In oiu* case 

T, (^)^"-"^-. 7'„(^)^"""'" = y^^«©'"""'"= const.; 

hence the theorem holds true and the critical value of 
Pm is determined from the condition 

(n-l)/n /^.\(n-l)/n 



\P2/ \Pml 



whence 

, /T \in/(n-l) 

P«. = Vpip2 • (|f j .... (23) 

If it is possible to cool the air to a temperature Tm = T2 
the preceding formula is simpUfied to 

Pm = Vpip^ (24) 

Prob. 16. — Show from Eq. (21) that when Pm satisfies Eq. (23) 
the work done in the low-pressure cylinder is greater than that done 
in the high-pressure cylinder when Tm > ^2, and vice versa. What 
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conditioli should Pm satisfy in order that equal amounts of work be 
done upon the gas in both cylinders? 

Prob. 17. — Prove result (23) by directly differentiating Eq. (22). 

Prob. 18. — Solve Prob. 15 assuming the machine to be a double- 
stage compressor, the air being cooled in the receiver to a temper- 
ature Tm = 40® C; determine the per cent saving in power. 

Ans, Pfn = 3.79 metric atmos. absolute; 

(a) Low-pressure cyl. 0.0925 cu. m.; high-pressure, 0.0262 cu. m. 
(6) " " 46.43 kw.; " " 35.52 kw. 

(c) 107** C. 

(d) " " 65.7% " " 54.0% 

Per cent saving = 8.44. 

Multi-stage Compressors. 

When large quantities of air are to be compressed 
to very high pressures, and economy in operation is of 
importance, the work of compression is divided among 
three or more cyUnders so that the air may be cooled 
in intermediate receivers. The problem is to deter- 
mine the pressures in these receivers so that the total 
work of compression may be a minimum. By analogy 
with Eq. (22) we have 

(^ \(n-l)/n /^ \(n-l)/n / ^r^^ \(n-l)/n 

^) +^4?") • • '+^A^) (25) 

P2 / APm/ \Pmk/ 

there being k intermediate receivers, or k + 1 compres- 
sion cylinders; the pressxu-es and the temperatures in 
these receivers are denoted by subscripts mi, ?n2 , . . . , 
m,. All the k mtermediate pressures are independent 
variables; all the temperatures are supposed to be 
given. According to the general rule of differential 
calculus the k partial derivatives of y with respect to 
the k independent variables must be equated to zero; 
but here again advantage may be taken of the theorem 
proved in Prob. 21 below, since the product of all the 
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terms constituting y is constant. Hence, all the terms 
must be equal to each other, or 

(^ \(n-l)/n /^ \(n-l)/n / 7), \ (»"!)/» 

From these k equations the k unknown pressures may 
be easily determined. 

Prob. 19. — Prove that when all the temperatures in Eq. (26) 
are equal, the pressure in each intermediate receiver is a geometric 
mean of the pressures in two adjacent receivers. In other words, the 
pressures increase according to a geometric progression whose ratio is 



Prob. 20. — To find a minimum value of the expression 

y = u + v, (27) 

where both u and v are functions of the same independent variable 
Xj and these functions are such that 

uv = const (28) 

Solvdon, The condition for a critical value of y is 

dy = du + dv = 0, 
while from Eq. (28) 

udv + vdu = 0. 

Eliminating du and dv from these two equations the result is 

u = v (29) 

It is left to the student to show that the values of u and v so 
foimd actually correspond to a minimum and not to any other 
critical value of y. 

Prob. 21. — Referring to the preceding problem, let, in a more 

general case, 

y = ui + U2 + uz+ - ' ' +Uk, . . . . (30) 

where each of the u*8 is a function of n independent variables xi, 
X2, . . . , Xn. Let also the condition be fulfilled that the product 

uiU2Uz , . . Uk = C = const. . . . (31) 
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It is required to find the values of the u*s for which y becomes a 
Tnininmim . SolvMon: Function y has a critical value for the values 
of xi, X2, . . . , Xn determined from the equations 

ia = 0;^=0;^=0;...-^ = 0. 

dXi dX2 dXz dXn 



In our case these conditions become 

1 1 1- . , . 

dXi dXi dXi 

dUi . dlh , dUz , 
dX2 dX2 dX2 



dxi ' 

dX2 



&Ui . dtt2 I ^m I 
dXn dXn dXn 



• • • 



dXn 



. . (32) 



The condition (31) being differentiated in succession with respect 
to the various x's, gives 



1 dl^i , 1 dll2 , 

— • 1 — • \- 

Ui dXi U2 dXi 
1 dl^i , 1 dU2 . 

— . 1 • 1- 

Ui dX2 U2 dX2 



Uk dXi 

+ 1.^ = 0, 

Uk dX2 



Ui dXn U2 dXn 






(33) 



The two systems of equations must be satisfied simultaneously. 
This is possible only when 

Ui = U2 = ui= • • • = wjb = C^/* (34) 



This is the required condition for min. y, although, strictly speak- 
ing, the values of y in the vicinity of the critical value must be 
investigated before one may state positively that condition (34) 
corresponds to a minimum and not to some other critical value of y; 
see Prob. 22. 

To prove that Eq. (31) being differentiated with respect to X\ 
gives the first Eq. (33): The partial derivative of Eq. (31) with 
respect to xi is 
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{UiUz . . . Uk) -^ + (UiUz . . . Uk) — 

dXi dXi 

+ • • • + (UiU2 . . . Uk-i) -7^ = 0, 

dXi 

or 

C dUi . C dU2 , , C dUk n 

Ui dXi Ui dXi Ujc dXi 

Cancelling C the first equation (33) is obtained. 

Prob. 22. — Prove that the values of w's according to Eq. (34) 
actually correspond to a minimum and not to a maximum of function 
y in Eq. (30). Solution: Solve Eq. (31) for Uk and substitute its 
value in Eq. (30). This gives 

jfc-i 
2/ = 2) ^i + C/(uiU2 ' • • Uk^i) .... (35) 
1 

AU the u^s in this equation are independent of one another, so 
that for the critical value of y all the first partial derivatives of y 
must be separately equal to zero. Taking the derivative with 
re8i)ect to a variable it<, we get 

^=1 - = 0, ... (36) 

dUi UiUi ' ' • Uk-1 Ui 

or 

uiiii • • • Uk-i • Ui = C (37) 

Comparing this expression with Eq. (31) we see that 

Ui = Uk (38) 

Since the subscript i has all the values from 1 to A; — 1, this result 
is identical with Eq. (34), and shows a simpler method of deriving 
it, by eliminating one of the variables. 

Differentiating Eq. (36) again, we find two kinds of second partial 
derivatives: 

d'y _ 2C _ 2 _ 2 ^ ,3^. 

dU^ U\Ui • • • Ujer-l • U^ Ui C^/* 

and 

-^ e ^ 1 = J- . . (40) 

dUi&Uj UiV/i ' ' • Ujt-l • U^j Ui C^/* 

Now let each variable receive some increment with respect to its 
critical value C*/*, and let this increment for w< be represented in 
the form C2»-, where e is the &ime for all the w's, but z is different 
for each u. If we can prove that the difference between the new 
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and the old value of y is positive for any small values of €«<, then it 
will follow that y is a minimuni for the critical values of the variables 
according to Eq. (34). Applying Taylor's series we have 

But at the critical point all the first derivatives according to Eq. 
(36) are equal to zero, and the values of the second derivatives are 
given in Eqs. (39) and (40). Consequently, 

^V = ht^W + ^2^+ • • • + 2jb-l2 + (21 + «2 

+ • • • + 2fr-i)*]/Ci/* + etc. . . . (42) 

The factor € can now be taken so small that the term containing e^ 
is larger than all the succeeding terms put together. The expres- 
sion in the brackets in Eq. (42) being a simi of squares is always 
positive; hence Ay is always positive whatever signs the Zi have. 
Therefore y is a TninimnTn when the variables are all equal to each 
other, according to Eq. (34).* 

* The author is indebted to Dr. Virgil Snyder, professor of mathe- 
matics at Cornell University, for this proof. 



CHAPTER III. 

PRESSURE-TEMPERATURE CURVES OF 
SATURATED STEAM. 

The action of the steam engine and. of the steam 
turbine is based upon the elastic properties of steam; 
therefore it is important for a research engineer to 
know how to represent these properties mathemati- 
cally. One of the most important data on saturated 
steam is the relationship between its temperature and 
the pressure at which it is formed from liquid water. 
The underlying phenomenon is 
as follows (Fig. 12): Let a unit 
weight of water be contained in 
a cylinder C, covered with a 
piston. The piston is loaded 
with a weight P such as to pro- 
duce upon the liquid a pressure 
of p kg. per sq. m. The water 
" is gradually heated until steam 
begins to form imder the piston. 
The temperature is then noted 
on the thermometer T. If more 
heat is added to the water, more 
am is formed, and the piston 
^ns to rise. By observing 
the thermometer it will be noted that the temperature 
now remains stationary until the last drop of water 
has been evaporated. After this, the temperature of 
the steam begins to rise, and it gradually acquires the 



Chap. III.] PRESSURE-TEMPERATURE CURVES. 41 

properties of a perfect gas considered in Chapters I and 
II. As long as liquid water is present the steam is 
called saturated and its temperature depends only upon the 
pressure p. Adding heat to it does not raise the tem- 
perature of the vapor, but merely increases the quan- 
tity of steam and decreases the amoimt of water. 

So long as any water is present, steam is called wet; 
when the last drop of water has been evaporated the 
steam becomes dry; if more heat is added the steam 
becomes superheated. 

The only way to obtain saturated steam at a higher 
temperature is to increase the pressure p. By repeat- 
ing the same experiment with a larger weight P it will 
be foxmd that the water must reach a higher temper- 
ature before steam begins to form. Thus, at ordinary 
atmospheric pressure, water begins to boil or steam is 
formed at 100° C. At a pressure of 5 (metric) atmos- 
pheres the temperature of water must reach about 151° 
C. before steam begins to be developed. 

The experimental relationship between the pressure 
of saturated steam and the corresponding temperatiu^e 
is shown on the curve sheet Fig. 13. The true theoret- 
ical equation of this curve is unknown, although quite 
a number of empirical formulse have been proposed 
which express this curve with a suflGicient acciu'acy. 
For ordinary practical pmposes no formulae are re- 
quired, values of pressures for different temperatiu-es or 
vice versa being taken from numerous " steam tables," 
which represent corrected and interpolated experimental 
results. But there is one important property of sat- 
urated steam, namely, its density, which cannot be 
very accurately determined experimentally but which 
can be easily calculated from the shape of the p-t curve. 
Moreover, in some theoretical investigations involving 
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differentiation and integration, it is necessary to have 
an analytical expression for the p-t curve of saturated 
steam. 

Thus, the problem in this chapter is to investigate 
various analytical expressions for the curve given in 
Fig. 13 and to show methods for calculating the nu- 
merical coefficients entering into such formulae. In 
doing so we shall limit ourselves to the range of temper- 
atures between 100 and 200° C, as the most important 
in engineering work. The range between 40 and 100° C. 
is also of importance in condensing engines and low- 
pressure turbines, but the above range is sufficient in 
order to illustrate the mathematical method.* 

In attempting to express a phenomenon by an em- 
pirical formula the simplest mathematical functions 
are tried first. Thus, there are two classes of functions 
or curves which have the general shape suggested by 
Fig. 13, namely paraboUc curves of the form 

p = Ax", (1) 

and exponential curves of the form 

p = 06*' (2) 

The first curve would mean in our case that the pressure 
p increases as the nth power of the temperature; the 
second curve would mdicate that the pressure increases 
in a geometric progression, when the temperature in- 
creases according to an arithmetic progression. Nu- 
merous empirical formulas have been proposed of the 
form (1) or (2), and some of them are considered below. 

* For an exhaustive treatment of the subject from the point of view 
of theoretical thermodynamics see G. A. Goodenough, ''Thermal 
Properties of Steam/' Bulletin No. 75 of the University of IlUnois, 
1914. 



44 ENGINEERING MATHEMATICS. [Chap. III. 

Parabolic Formulae. 

First assume the curve in Fig. 13 to be an ordinary 
parabola of the second degree with its vertex some- 
where to the left of the curve sheet. Let the unknown 
pressure at the vertex be po and the temperature to. 
With respect to this point as the origin, the ordinates 
increase as the squares of abscissae, so that 

(p - Po) = A (^ - toy, .... (3) 

where A is a coefficient of proportionality, to be de- 
termined from the given curve. 

Eq. (3) contains three unknown quantities, to, po, 
and A. These may be determined by applying the 
equation to some three points on the given curve, for 
example, to the points corresponding to the temperatiu'es 
of 100, 150 and 200° C. Denoting these temperatures 
fi, fe, and <3, and the corresponding pressures by pi, p2, 
and Pa, and assuming that the theoretical curve (3) 
passes through these three points on the experimental 
curve, we have 

pi -po = A {ti - toy, .... (4) 

P2 - Po = A (<2 - toY, .... (5) 

P3 - Po = A (fe - toy (6) 

These three simultaneous equations may be solved 
with respect to po, to, and A. By subtracting Eq. (4) 
from Eq. (5), and Eq. (4) from Eq. (6) two new equa- 
tions are obtained which do not contain po. Dividing 
one of these equations by the other the unknown quan- 
tity A is eliminated, and finally to may be determined. 
Expression (3) is of the first power with respect to 
p and of the second power with respect to t It is 
therefore equivalent to the equation 

p = a + bt + cP, (7) 
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in which the constants a, b, and c can be determined 
by writing Eq. (7) for some three points on the experi- 
mental curve (Fig. 13) and solving the three resulting 
simultaneous equations. 

The actual shape of curve (3) or (7) depends on the 
three points selected on the given experimental curve. 
The best theoretical curve is not necessarily the one 
which crosses the experimental curve at the points 
corresponding to the temperatures of lOO"^, 150°, and 
200° C. The best theoretical curve, or, as it is usu- 
ally called, the most probable ciu^ve, may cross the 
experimental curve in some other points. Therefore, 
the above method does not give the best possible solu- 
tion for the selected general form of Eq. (3) or (7). 
One way to improve the numerical coefficients is by 
inspection and trials; a more accurate procedure is to 
apply the method of least squares which is explained at 
the end of the chapter. 

Prob. I. — Obtain expressions for po, fc, and A as explained above, 
plot curve (3), and determine the maximum per cent error or devi- 
ation from the experimental curve in Fig. 13. 
Ans. With the data corresponding to < = 100, 150, and 
2()0° C. the values are: po = 1.0263 kg. per sq. cm.; 

«o = 98.4° C; A =0.001442. 
Prob. 2. — Express the constants a, 6, and c through'po, ^o, and A, 

Ans, a = (po + -^V); h = — 2 A^; c = A, 
Prob 3. — Find the numerical values of the constants, a, 6, and 
c, by solving the three equations which represent relationship (7) 
for some three points on the experimental curve. 
Ans, For t = 100, 150, and 200° C, a = 15; 6 = -0.2839; 

c = 0.001442. 
Prob. 4. — Solve the preceding problem in the general form, 
that is, find algebraic expressions for a, 6, and c, when the parabola 
(7) is to pass through three given points on the experimental curve, 
these points being characterized by abscissae ti, ^, and tsj and 
ordinates pi, p2, and ps. 
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The next parabolic formula to try is one of the form 

p = A<~ (7a) 

This formula can be made to hold true only for a 
limited range of the curve, say between 100 and 200° C. 
It differs from Eq. (3) in that it represents a parabola 
of the nth degree, where n is to be determined from the 
experimental curve. Moreover, the constants po and to 
are at first assumed equal to in order to make the for- 
mula simpler. Taking common logarithms we get 

log p = logA +nlogt. . . . (7b) 

Thus, plotted to a logarithmic scale for p and t, the 
curve becomes a straight line of the form y = a + nx. 
This fact suggests a simple method for determining A 
and n graphically: plot the given curve (Fig. 13) on 
a sheet of logarithmic paper and draw a straight line 
which comes nearest to the ciu^e. The trigonometric 
tangent * of the angle which the curve makes with the 
axis of abscissae gives the exponent n. Knowing n, the 
parameter A is determined from Eq. (^b) by substitut- 
ing a pair of known values for p and t 

A paraboUc formula which fits the experimental 
curve more closely is 

p = A(t-toy (8) 

It differs from Eq. (7a) in that it contains an additional 
constant k- The coeflficient po is again assimied equal 
to 0, in order to keep the niunber of constants down to 
three. A more rational and general formula would be 

(p - po) = A (< - to)% .... (9) 

* The trigonometric tangent represents the actual slope of the line, 
and should be measured by an ordinary scale and not by the logarith- 
mic scale. 
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which contains four coeflBicients and therefore can be 
made to pass through any four desired points on the 
given experimental curve. Applying Eq. (8) to some 
three given experimental sets of data we get : 

Pi = A {ti - uy, 

P2 = A(<2-fe)«, . . . \ . (10) 
pz = A(k- <o)^ 

From these equations A is easily eliminated by division, 
the result being 

2>2 \fe — to/ ' 



2? ^ / fe - to Y 

Pi \ti - to) ' 



(11) 



From these two equations n may be eliminated by tak- 
ing logarithms of both sides of the equations and divid- 
ing one result by the other. The resultant equation for 
to is 



log 7>3 - log P2 _ 
log p2 - log pi 






. . (12) 



This equation can be solved by trial only, assuming 
different values for to, beginning, say, with zero. The 
procedm-e is not tedious, since the left-hand term of 
the equation is a known quantity. If in Eq. (12) the 
value of p2 is so selected that 



(13) 



23 ^ 2? 

P2 Pi' 

the left-hand side of the equation becomes equal to 
unity, and consequently 

tz — ^0 fe — to 



ti — Co fe — to 



(13o) 
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From this equation to may be found. Or the values 
of pi, p2, Pz may be so selected that 

^ (14) 



©- 



P2/ Pi 

This makes ttie left-hand side of the equation equal to 
2f so that a quadratic equation for to is obtained. 

Instead of eliminating n, the coefficient to may be 
eliminated from Eq. (11), by extracting a root of the nth 
power and determining to from both equations. The 
resultant equation for n must be solved by successive 
approximations imless the value of p2 can be so selected 
as to satisfy Eq. (13) or (14), in which case n may be 
computed directly. 

Eq. (8) plotted to a logarithmic scale gives a straight 

line 

logp = log A + nlog {t - to), 

between p and {t —to). Therefore, in order to see 
how closely a paraboUc curve checks with the given 
experimental curve, it is sometimes convenient to 
plot both curves on a sheet of logarithmic paper for 
comparison. 

Prob. 5. — Plot the curve in Fig. 13 to a logarithmic scale and 
determine A and n in Eq. (7a) . Check a few points with the original 
curve. Am. A = 0.127 X 10~^; n = 3.95 (approx.). 

Prob. 6. — Determine the numerical values of the coefficients 
in Eq. (8) first using Eqs. (12) and (13), then Eqs. (12) and (14). 
Ans. For pi = 1.03 and pz = 15.9 kg. per sq. cm., by Eqs. 
(12) and (13), to = -48.6° C, by Eqs. (12) and (14), 

to = -39.1*' C. 
Prob. 7, — Determine the numerical values of the coefficients 
in Eq. (8) so that the curve would pass through the points corre- 
sponding to ti = 100**, t2 = 150°, tz = 200° C. HirU. — Solve for 
to first by trial. 

Ana, to = -39° C.]n = 5.06; A = 0.149 X 10-" 
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Prob. 8. — Eliminate ^o from Eqs. (11) and show how to deter- 
mine n, A, and to. Also simplify the method for the case when 
either Eq. (13) or (14) is satisfied. 

Prob. 9. — Outline the method by which the coefficients in for- 
mula (9) may be calculated. Determine roughly their numerical 
values, and check the results on a sheet of logarithmic paper. 

If formula (7) does not represent experimental results 
with sufficient accuracy, it may be improved by adding 
further correction terms dP + et^ +ft^ + ' ' ' . Or 
the terms added may be of the form 

Trial alone can tell which form gives a better approach 
to the given ciu*ve. If it is desired to limit the nmn- 
ber of terms to three, fractional exponents may be used, 
so that the formula becomes 

p = a + bV^ + ct\ .... (15) 

This formula contains five constants, and can there- 
fore be made to pass through any five points on the 
given experimental curve. It can be made to fit it 
more closely than a curve with only three coefficients. 

Another convenient formula with five coefficients is 
of the type 

"-^Irf <•« 

It has the advantage over Eq. (15) in that the pressure 
p may be more conveniently calculated by means of a 
table of logarithms. For vapors of some substances 
satisfactory results are obtained by using only four 
instead of five constants in Eq. (16) ; that is, one may 
assiune a = 0, or j8 = 0, or m = n. 
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Prob. 10. — Show how to calculate the coefficients o, 6, and m 
in formula (15), aBsuming c = 0. 

Prob. zz. — Same as Prob. 10, only c is not equal to zero. 

Prob. Z2. — Show how to determine the coefficients in formula 
(16) when m = n. Hint, — Select on the experim ental curve such 
points^t pi : p2 = P2:pz = Pi :p4, and denote Vp^ = </pz/P2 
= ^/pi/pi by X ; this quantity can then be eliminated. The problem 
is thus reduced to a solution of two simultaneous equations for a 
and/9. 

Logarithmic Fonziule. 

We now turn our attention to the other class of 
formulae of which the simplest is Eq. (2). Taking 
common logarithms the formula becomes 

logp = A+Bty (17) 

where A and B are constant coeflScients. If this for- 
mula expressed correctly the relationship between the 
pressure and the temperature of saturated steam, the 
experimental curve (Fig. 13) should become a straight 
line when plotted to a logarithmic scale for p and to 
a uniform scale for t. This is only approximately true, 
but may be accurate enough for some purposes. 

Formula (17) may be improved by adding terms to it; 
for instance, 

logp^A+Bt + Ct^+ (18) 

Other self-evident modifications are 

Bt 

logp=A+^^; .... (19) 

logp=A-|-|; .... (20) 

logp = a+^+7log<, . . . (21) 

etc., etc. 
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Formula (19) was used to a considerable extent in 
the early days of steam engineering. Formula (20) 
is known as Rankine's formula. As a rule, absolute 
temperatures are used in it. 

Pjob. 13* — Replot the curve in Fig. 13, using a logarithmic 
scale for ordinates; draw a straight line that comes closest to the 
curve, and determine the constants A and B in formula (17). 

Prob. 14. — Determine the numerical coefficients in formula (20), 
using absolute temperatures for t. 

Am. A = 7.63; B = 3940; C = -423,600 for 

T = 373° C, 423° C, and 473° C. 

Formula of Biot and Regnault. 

Regnault, to whom we are mainly indebted for our 
first knowledge of the properties of saturated steam, 
was confronted with the problem of presenting the 
relationship in Fig. 13 by a suitable formula. The 
simplicity of the formula was of secondary consider- 
ation, extreme accuracy being of prime importance. 
He selected the form known as Biot's formula 

log p = A + Bb' + Cc\ . . . (22) 

which is somewhat similar to Eq. (18), except that the 
temperature enters as an exponent. Regnault probably 
found that this form of equation satisfied the results, of 
his experiments better than any other forms known to 
him.* The curve corresponding to Eq. (22) may be 
made to pass through five points on the experimental 
curve, because the equation contains five coefficients; 
that is, Eq. (22) may be written for five known sets 

* The original form given by Regnault was 

log p = A + ik6*-*o 4- fc*-<», 

where ^ is the lowest limit of the range of temperature to which the 
formula applies. Denoting kh~^ by B and hr^^ by C expression (22) is 
obtained. 
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of values of p and t, and the five simultaneous equations 
thus obtained may be solved for A, ByC, 6, and c. 

In the general caee these five equations can be solved 
by trial only, which is an extremely tedious procedure. 
However, an exact solution is possible if the five pres- 
siu'es p are selected to correspond to equidistant tem- 
peratiu*es; in other words, when 

14 — tz =^ tz — fe = fe — ti^^^i — M) = 5. . (23) 

The method for determining the constants in this case 
is as follows: Applying Eq. (22) to five equidistant 
points on the curve, we get 

logpo = A + Bbh +Cch, 

Jog pi = A + B6^+« +C(^-^', 

log p2 = A + By«+2a + c^+2a^ [ ... (24) 

log p3 = A + Bi/'+^j _^ Cc^+3«^ 

log p4 = A + Bb^-^' + C&^-^\ 

These five simultaneous equations must be solved with 
respect to the five unknown coefficients Ay Byb, C, c. 

Subtract the first equation from the second, the 
second from the third, etc., in order to eliminate A. 
The result is 

log pi - log po = BV^ (6* - 1) + Cc^ (c* - 1), 
log p2 - log pi = Bb^^' (6* - 1) + Ccf^-^' (d - 1), 

log pz - log P2 = Bb^+2a (^a _ 1) + Cc^+2a (c« - 1), 

log p4 - log Pz = BV^+^' {¥ - 1) + Cc^+3*(c* - 1). 

Now denote 

B6'«(6f- 1) =x, 

C(f^{c''-1) =y, 
V = V, (f = w, 
log pi - log Po = go, [ . . . . (25) 
log p2 - log pi = gi, 
log Pz - log p2 = g2, 
log p4 - log p3 = qz. 



r" 



(26) 
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Then the four preceding equations assume the simple 

form 

qo = x + y, 

qi = XV + yw, 

32 = xv^ + yw^, 
qs = xt^ + yw^. 

The standard method for solving such equations is to 
determine x from the first of them and to substitute 
its expression into the three remaining equations; 
then to determine y from one of the equations so ob- 
tained and to substitute its value into the remaining 
two. The result is two equations containing the un- 
known quantities v and w, which can easily be de- 
termined. 

There is, however, a shorter way to eUminate x and 
y simultaneously, namely, the so-called method of 
indeterminate coefficients. Let us assume that one 
of the equations between v and w is the result of an 
elimination of x and y from the first three Eqs. (26), 
while the second equation between v and w is the re- 
sult of the elimination of x and y from the last three 
Eqs. (26). Multiply the first Eq. (26) with an inde- 
terminate coefficient P, the second with a similar 
coefficient Q, and add them to the third equation. 
The result is: 

qoP +qiQ +q2=x{P+vQ +v^) +y (P +wQ +w^) . (27) 

This equation is true for any values whatever of P and 
Q. As our purpose is to eliminate x and y, we shall 
give P and Q such values that the coefficients of x and 
y shall become zero; in other words 

P + vQ +v^ =0, 
P + wQ+w^ = 0, 



54 ENGINEERING MATHEMATICS, [Chap. IH. 

from which 

qZ'^Iv+w)} ^^^ 

Thus, Eq. (27) becomes 

Qo ' vw — Qi {v +w) + q2 = 0. . . (29) 

It is not necessary to repeat this process for the 
last three Eqs. (26), because the equation obtained by 
their addition is identical with Eq. (27) except that xv 
stands for x and yw for y. Therefore, the values of P 
and Q are the same, and we can write directly 

Qivw — q2{v + w) +qz = 0. . . . (30) 

Eqs. (29) and (30) can easily be solved for v and w. 
The simplest way is to consider vw to be one unknown 
quantity and {v + w) to be the other. Using the 
notation of Eqs. (28) we have 

qiP + Qq2+qz = 0.y * * ' ^^ ^ 

These must be solved for P and Q, and then Eqs. (28) 
solved for v and w. Knowing v and w, the unknown 
coefficients 6, c, B, C are determined from Eqs. (25), 
and the problem is thus solved. 

Prob. 15. — Show that by a direct elimination of x and y from 
Eqs. (26) results are obtained identical with (29) and (30). 
Prob. 16. — Solve Eqs. (31) for P and Q. 

Ans, P = (qtqi - 52*)/(Mo - ?!*), 
Q = (Mi - ««9'o)/(Mo - gi*). 
Prob. 17. — Solve Eqs. ( 28) for v a nd w. 

Am, w;=i(-Q±^Q'-4P); t;=i(-Q=f=^Q'-4P). 
Prob. 18. — Determine the coefficients in fonpula (22) from five 
equidistant points on the curve in Pig. 13. Use five-place loga- 
rithms and be careful about the last significant figures.* Note, — 
5 must be odd in order to arrive at a real value for b or c, 

* Regnault used Vega's ten-place table of logarithms in the calcula- 
tion of these coefficients. 
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Prob. 19. — Regnault found that at low temperatures and pres- 
sures the last term in Eq. (22) may be neglected. Show how to 
determine the coefficients A, B, and b in this case. 

The Most Probable Curve. 

With the method outlined above, the exact values of 
the coefl&cients in a formula depend upon the selection 
of the points on the experimental curve through which 
the approximate curve is supposed to pass. Thus, 




1— Experimental 

2 — Passing through ends 

8-.- Most probable 



Fig. 14. — The most probable curve as distinguished from a curve 
which passes through the end points of the experimental curve. 

assuming, for instance, the form of the equation to be 
that represented by formula (7) the coefficients a, 6, 
and c would probably be determined so that the curve 
would pass through the extreme and the middle points 
on the given cm^e (curve 2, Fig. 14). But this may 
not be the best solution. It is possible to modify 
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sligkUy the coeflBicients a, 6, and c in sucH a way as to 
make Eq. (7) check much more closely with the given 
curve (see dotted curve 3 in .Fig. 14). 

Upon first thought, the most probable curve seems 
to be the one for which the sum of the errors is a 
minimum. An " error " is here defined as the differ- 
ence of the two ordinates corresponding to the same 
abscissa, one ordinate belonging to the experimental 
curve and the other to the assumed or calculated curve. 
But a closer reasoning will show that the sima of the 
errors for curve 2 (Fig. 14) may not be larger or may 
even be smaller than the sum of the errors for curve 3, 
because positive errors are balanced by negative errors. 
And yet it. is evident that curve 3 satisfies the experi- 
mental results better than does curve 2. 

It is shown in the theory of probabiUties that for 
the most probable curve the sum of the squares of the errors 
is a minimum. One can readily see that whether an 
error is positive or negative its square is always positive ; 
if the absolute values of errors are reduced by making 
the curve to pass more closely through the experimental 
points, the squares of the errors are reduced to a still 
greater degree. It is reasonable to assmne that the 
curve for which the smn of the squares of the errors 
is a minimum is the best possible curve. 
. The above law of least squares is firmly established 
in application and one of the arguments in its favor 
is that in the simplest case of one quantity to be de- 
termined, the method leads to the arithmetical mean 
of the observations. For instance, let the distance 
between two fixed points be measiu-ed k times, in order 
to obtain a greater precision than is possible with a 
single measurement. Let the results be Zi, tj, fe, • • ., 
Zjfc, which differ slightly from one another. Let I be 
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the most probable unknown value of the required dis- 
tance. Then, if the law of the least squares is correct, 
the expression 

S (errors)^ = (I - ky + {I - k)^ + • • - + {I - hY 

must be a minimum. Taking the first derivative of 
this expression with respect to the imknown quantity Z, 
and equating it to zero, we have 

2 (Z - ii) + 2 (Z - ^2) + • • • + 2 (Z - Z*) = 0, 
or 



I = 



Zl + ^2 + • • • Zfe 

k 



In other words the most probable value of the distance 
is the arithmetical mean of the observations, as is sug- 
gested by conunon sense. Perhaps this proof also 
gives an answer to the question as to why a sum of the 
squares is taken and not a sum of any other even powers 
of errors. 

The simplest case of application of the method of 
least squares is found when the dependent variable, in 
our case p, and the coeflScients to be determined, such as 
a, 6, and c, in Eq. (7), enter in a linear form. Take any 
number of points on the experimental curve, say, twenty 
or thirty, and let them be determined by the abscissae 
^1, fe, . . . , k, and by the ordinates Pi, P2, . . . , p*. 
The theoretical curve (7), as a rule, does not pass 
through any of these points, so that Eq. (7) becomes 



Pi + Ai = a + 6<i -h cti\ 
P2 + A2 = a + 6fe + cti^, 



p* + A* = a + &<A + cti?, 



(32) 
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where Ai, A2, . . . ; Aib are small errors between the 
values of p calculated from Eq. (7) and those actually 
observed. 

Accordmg to the law of least squares, if the coeflS- 
cients m Eq. (7) have been properly selected so as 
to give the most probable curve, the errors satisfy the 
condition 

Ai2 + A22 + • • • + Afc2 = min., 

^^ X^^ + ^^i + ^*^ - P*)' = ^i^- • • • (33) 

In this equation all the p's and^'s are taken from the 
experimental curve. The only imknown variable 
quantities 'are the^coefficients a, 6, and c. Since these 
coeflScients are independent of each other, function (33) 
becomes a minimum when each of the three partial 
derivatives 

dS dS as 

da' db^ dc 

is separately equal to zero. 
Performing ^the differentiation and dividing by 2, 

we get 

2 {a + bt + cP - p) = 0, 1 

2< (a + &< + d2 - p) = 0, . . . (34) 
2<2(a + W + c<2-p) =0, 

where the subscripts i are omitted for the sake of 
brevity. Simplifying, we obtain: 

ka + 6S< + cSP = Sp, 
dZt + &2<2 + cS<« = 2p«, ... (35) 

Solving these simultaneous equations for a, 6, and c, the 
most probable curve is obtained. It will be observed 
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that when usmg the method of least squares it is not 
necessary to specify in advance through which points 
the curve should pass, nor are we limited in the num- 
ber of experimental points to be used in the solution 
of Eqs. (36). 

r 

Prob. 20. — Determine the coefficients A and B in Prob. 13, by 
the method of least squares, and see how close the curve thus de- 
termined comes to the one drawn according to the judgment of the 

eye. 

Ana, For the ten points, t = 100, 110, 120, . . . , 190** C, 

A = -1.174,5 = 0.1224. 

A Method for Improving Empirical Formulae. 

A weak point in the above-described appUcation of 
the method of least squares is that the coefficients 
which enter into equations such as Eq. (35) are usually 
large numbers. Thus, in our case t varies between 
100 and 200, so that t^ is of an order of magnitude of 
thousands of millions. DeaUng with such large num- 
bers is not only inconvenient, but may cause inac- 
curacy in multiplication and division. Thus, an error 
of one per cent in any of the 2's introduces a corre- 
sponding error in the magnitudes of a, 6, and c. At 
the same time, unless a, &, and c can be determined 
with considerable accuracy, the whole method of least 
squares is superfluous, since approximate values of 
coefficients can be determined by simple trials. There- 
fore, it is preferable first to calculate some approxi- 
mate values of the coefficients in question without the 
use of the method of least squares, and then to apply 
the method of least squares to finding corrections to 
these approximate values which will make them more 
accurate. This procedure is technically known as 
" improving '' empirical formulae. An example will 
make it clearer. 
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Let the best numerical value of a coeflEicient in an 
empirical formula be 28.45(7), the digit in the paren- 
theses being uncertain. If this value be calculated 
by the method of least squares, from equations such 
as (34), computations must be conducted with an 
accuracy of at least six significant figures, the last one 
or two being lost during the transformations. Let, 
however, the approximate value of the same coefficient, 
say, 28.1, be determined by some other method, for 
mstance, by assigmng three points through which 
the theoretical cm^e is supposed to pass. Then the 
method of least squares need be appUed only to the 
correction 0.35(7) of this value (in the manner indicated 
below) ; it is sufficient to compute with, say, only four 
significant figures. Moreover, with such a procedure, 
some of the coefficients in Eq. (35) assume much more 
convenient values and are easier to handle. 

Let A, B, and C be some known approximate values 
of a, &, and c, and let a, j3, and 7 be small unknown cor- 
rections of these values, so that 



a = A + a, 
& = B + ^, 

c = C + 7. 



(36) 



Substitute these expressions in the first Eq. (32); 
the result is 

pi + Ai = A + a + B^i + iSfi + W + 7<l^ 
or 

ri+Ai = a + /3fi + 7<A (37) 

where 

ri = pi - (A + Bti + Ck^). . . . (38) 

Here n is a small quantity, being the difference between 
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the observed pi and that calculated from the approxi- 
mate formula p = A + Bt + Ct\ 

Thus, Eq. (37) is identical in form with Eq. (32), 
except that small quantities r, a, j8, and y enter in place 
of large quantities p, a, 6, and c. By identical reason- 
ing with that given above, we arrive at the conclusion 
that the corrections a, jS, and y must satisfy equations 
identical in form with Eqs. (35), namely 

ka + fi^t + yJlt^ = Sr, 

aUt + m^ + 72<3 = 2rt, . . . . (39) 

m 

Solving these for a, j3, and 7, and substituting the values 
into Eqs. (36), the most probable values of a, 6, and c 
are obtained. 

Prob. 21. — Determine approximate values of a, 6, and c in 
Eq. (7) so that the curve will pass through the points correspond- 
ing to i = 100, 150, and 200° C. Then calculate corrections accord- 
ing to the method of leaat squares [ase Eqs. (39)], selecting from 
ten to fifteen equidistant points on the curve. Plot the most 
probable curve on the same sheet with the experimental curve and 
with the curve of first approximation. 

Arts, a = 19.4; fi = -0.267; y = 0.00085; New equation 

therefore reads p = 34.4 - 0.551 1 + 0.00229 0, 

Non-linear Equations. 

The method outlined above cannot be used directly 
when the given equation is not linear with respect to 
the unknown coefficients. Take, for instance, Eq. 
(7a) : proceeding as in Eq. (32) we get 

Pi + Ai = Ati"", 

P2 + A2 = Afe**, 



p* + A* = Atk". 



(40) 
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According to the law of least squares the most prob- 
able values of A and n satisfy the equation 

S (Ar - p)2 = min. 

The partial derivative of this expression with respect 
to A being equated to zero gives 

sr (Ar - p) = 0. 

The partial derivative with respect to n gives j 

SArLnf (Ar - p) = 0, 

or, after reduction, 

AJlt^Lnt - ZpCLnt = J ' ' 

The unknown quantities A and n enter into these equa- 
tions in such a form that it is impossible to solve the 
equations. Such a case invariably arises when the 
coefficients do not appear in the given equation, such as 
Eq. (7a), in the linear form. The following artifice is 
used in this case in order to obtain linear equations in 
place of Eqs. (41). The coefficients A and n are again 
determined first approximately, for instance, as in 
Prob. 5. Let A' and n' be the approximate values, 
so that the most probable values are 

^^^Z + '^'j (42) 

where a and v are small quantities. Let the equation 
of the approximate curve be 

p' = A't-\ (43) 

while the equation of the most probable curve is 

p = At- = (A' + a) r'+' . . • (43a) 
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According to Taylor's theorem for two variables we 
have 

But, according to our supposition a and v are small 
quantities; therefore all the terms of the series, with 
the exception of the first three, may be neglected. 
Substituting the values of p' arid its derivatives from 
Eq. (43) we have approocimately 

Hence, the first of the Eqs. (40) becomes 

pi + Ai = A V + ctti^' + V . A V'Lnfi, . (45) 

or, introducing 

n = pi - A'<i»', 



we obtain 



Ai = at{'' + yA V'Ln^i - ri, 
A2 = ah''' + j'A V-Lnfe — T2y 



^k = cLtk''\ +yA VLn<* - n. 



(46) 



These are Unear equations with respect to the un- 
known corrections a and v. Therefore, applymg the 
condition SA^ = min., and taking partial derivatives 
with respect to a and v, we obtain two linear equations 
from which a and v are easily calculated. 

To the beginner this method may seem somewhat 
circuitous, but it is based upon the simple theorem 
that when the coefficients in a fimction receive small 
increments, the increase in the fimction itself is a 
linear function of these increments. In our case we 
change from the known function p', Eq. (43), to the 
unknown function p, by giving the coejfficients A' and 
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n' unknown increments a and v\ Eq. (44) indicates 
that the difference p — p' is a linear function of these 
increments. The difficulty with Eqs. (40) is that the 
right-hand side is not a linear function of the imknown 
coefficients A and n. By the abov6 artifice it is re- 
duced to the form (45) in which the right-hand side 
consists of a known function A'ti^' which need not be 
linear, and of a linear expression containing the im- 
known corrections a and v. 

Prob. 22. — From Eqs. (46) obtain two expressions for a and 
Vf which take the place of Eqs. (41). 

Ans, (1) a2<»«' + A'vi:^^'lM = 2r<~'; 

(2) a2<«»'Ln/ + AV2^«'(Lni)* = 2r««Xni. 
Prob. 23. — Using the results of Prob. 5 as the first approxi- 
mation, determine more accurate values of A and n by the method 
of least squares. Hint, — Use 4-place tables of logarithms for 
computing r, and be careful of the last figure. After obtaining 
two simultaneous equations for a and Vy use long-hand calculation 
for all changes, as the two equations differ very Uttle from each 
other, the first difference appearing in the fourth figure. 

Ana. a = 0.011 X 10-^; v = -0.018. 

New equation reads p = 0.1383 X 10"^ ^'^. 
Prob. 24. — Show how to improve the coeflficients in Rankine's 
formula (20) by the method of least squares. Solution. — The usual 
way would be to represent the formula in the form 

p' = 10 ' ^ 

and to expand it according to Eq. (44). - After this has been done 
formulse are obtained corresponding to Eqs. (45) and (46). The 
student is advised to follow this method and to obtain final expres- 
sions for the corrections a, /S, and 7, in order that he may be sure 
that he understands the general method. The same results can, 
however, be obtained in a somewhat simpler way by the following 
artifice: 

logp-log[p' + (p-p')l = logp'(l + ^') 



> > 
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where M = log c = 0.43429 is the modulus of common logarithms, 
and the expression in the last brackets represents the expansion of 
a logarithmic function into an infinite series. The difference be- 
tween the accurate value p and the approximate value p' is small, 
so that all the terms of the series may be neglected, except the 
first. Solving for p — p' gives 

p-p' = ^(logp-logp'), 
or taking into account expression (20), 

^ ^ M\ t ^1 
Hence, the expressions corresponding to Eqs. (40) become 



Pi+A. = p/ + ?i.«-^^--2L 



M 



Mti Mti 



2 '* 



P2+^2= P2 +'^a— ,*T7 /3 



22.^_ ^.__M. 



M 



Mti 



Mk^ 



(47) 



These expressions are linear with respect to the unknown correc- 
tions a, /9, and y. Denoting as before the difference pi — p' by 
fi, the law of the least squares gives: 

or, taking partial derivatives. 
After simplification we get 






t 



i'2 



-X^-^'X^-S^-^SI-^o- 



(48) 
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From these three equatioiis the ccnrections a, fi and y may be cal- 
culated. 

Piob. 35. — Ec}B. (48), from idiich the oofrecticMis aie to be cal- 
culated, are obtained in a definite manner from Eqs. (47). These 
latter equadoiis are also linear with leep&cX to the corrections a, 
fif jf etc., of the gieneral form 

Ai = aZ, + ^Fi + 7^1 + • • • - R, 

At^aXt-hfiYi + yZt+ • • • -Bt, 



Ajfc = aXk +fiYt + yZk+ • • • — Ri^ 

Deduce the rule according to which the coefficients in Ecjb. (48) 
are f onned from X,Y, Z, . , , , R, With such a rule there may 
be no need of actually differentiating eveiy time, in practical appli- 
cations, but equations such as (48) may be written down as socm 
as the problem is reduced to the linear form of Eq. (47). 



CHAPTER IV. 

HEAT OF VAPORIZATION AND DENSITY OF 

SATURATED STEAM. 

The phenomenon of vaporization of water under 
pressure, described at the beginning of the previous 
chapter, involves two quantities which are of prime im- 
portance in the design of steam engines and turbines, 
namely, heat of vaporization and specific weight or 
density of steam. Let us assmne that at the begin- 
ning of the experiment there is no steam in the vessel, 
and that the water is brought up to the temperature 
at which steam just begins to be formed (this tempera- 
ture depends upon the pressure exerted upon the piston, 
see Fig. 13). As soon as vaporization begins, the tem- 
perature remains constant until the last drop of water 
has been converted into steam. Hence, all the heat 
communicated to the mixture of water and steam is 
used for two purposes : 

(a) to raise the weight of the piston (external 
work) ; 

(6) to separate particles of the water in convert- 
ing it into steam (internal work). 

Experiment shows that at usual pressures and temper- 
atures the internal work is several times as great as 
the external. 

It will be remembered that in a perfect gas the dis- 
gregation work is zero, so that in an isothermal process, 
as is the one under consideration, all the heat received 

67 
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is converted into external work. In the case of a 
saturated vapor, however, many times more heat must 
be communicated than the amount converted into 
work. This heat remains in steam in a '^ latent state," 
and is Uberated again when the steam is condensed 
into water. Thus, 

Total heat of vaporization = A (external work) 

+ (internal latent heat), 

where A is the coefficient for the conversion of work 
units into heat units (see p. 9). The external work 
can be easily calculated since the expansion is per- 
formed at a constant pressure. 

External work = / pdv = p(v2 — vi), 

where Vi is the specific volume of water, t^ is that of 
steam. It is customary to denote (t^ — vi) by Uy or 

U = V2 — Vi, 

where u stands for the increase in the specific volume of 

saturated steam over that of water. For instance, 

at 150° C. one kilogram of saturated steam occupies 

a volume i^ = 0.3876 cu. m., while one kilogram of 

water at the same temperature occupies less than 

0.001 cu. m. Thus, at 150° C. u = 0.3876 - 0.001 or 

0.3866 cu. m. per kg. It will be seen from this example 

that for most practical piu'poses Vi may be neglected 

and u considered equal to t^. The above equation then 

becomes 

r = Apu + p, (1) 

where r is the total heat of vaporization, p is the in- 
ternal heat of vaporization, and the difference of the 
two, Apu, represents the external work done ; r, p, and 
u are, of course, values for unit weight of working 
substance. 
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Prob. I. — At a temperature of 150° C, steam has a pressure 
of 4.869 kg. per sq. cm., and r = 500.79 kg-cal. What per 
cent of total heat r is converted into work, and how much is re- 
tained by the steam as internal latent heat p? Hird, vi and vt 
are given above. Ans. 8.81 per cent; 91.19 per cent. 

Relation between r and u. 

Total heat of vaporization r, and specific volume u, 
can be determined separately by experiment. There 
is, however, a general theoretical relation between the 
two, so that when one of them is known the other may 
be calculated. This relation, known as Clapeyron's 
equation, may be deduced as follows: Let a imit 




Volume u 
Fig. 15. — A thermodynamic cycle performed by a mixture of 

water and steam. 

weight of water perform a thermal cycle indicated in 
Fig. 15. At A the water is all Uquid; it occupies a 
certain volume Vi, has a pressure p, and a temperatiu'e 
T (absolute) such that the water is about to vaporize. 
Heat is now being added to the water, which vapor- 
izes until at B we have dry saturated steam. Dur- 
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ing this process both the pressure and temperature 
remain constant. At B the pressure is reUeved by an 
infinitesimal amomit dp, and the steam aUowed to 
expand adiabaticaUy (see p. 15) along the line BC, its 
temperature dropping by dT. After this, the two 
processes are repeated in the reverse order; that is, 
first, heat is withdrawn at the temperature T — dT, 
so that the steam condenses into water along the iso- 
thermal CD, and secondly, the pressure is again in- 
creased by the amount dp, so that the substance is 
returned to its original condition along the adiabatic 
DA. Since the substance has now the same" heat con- 
tent as at the beginning of the cycle, the net amoimt 
of heat expended during the process is equal to the 
mechanical work performed during the cycle. The 
mechanical work is expressed on the diagram by the 
area ABCD, or is equal to 

{v2 — Vi) dp = udp (2) 

The total heat conmiunicated to the mixture during 
the isothermal expansion AB is equal to the heat of 
vaporization r. Of this only a small part is converted 
into work, the rest, r', being rejected at the temperature 
T — dT during the isothermal compression CD. Heats 
received and rejected under such conditions are as 
the corresponding absolute temperatures (second law of 
Thermodynamics); hence 

r T 
Qj. V'^ T -dT ^^^ 

^-f^'-f (4) 

Thus, the net heat communicated to the mixture is 

, dT 
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Equating this quantity of heat to the mechanical 
work performed, given by Eq. (2), we get Clapeyron's 
equation 

r— = At^dp (5) 

If r is the quantity to be determined, the equation is 
written in the form 

7=AuT-^, ..... (6) 
or if 16 is to be calculated 

^^Af'dp/dT ^^^ 

The statement contained in Eq. (4) may require 
some fiu-ther explanation. Let the dry steam, at 
point Bf be allowed to expand along Bp adiabatically 
(that is without addition or abstraction of heat) 
until the volimie becomes infinite, and the absolute 
pressure becomes zero. At that point, where all the 
adiabatics come asymptotically together, let compres- 
sion begin along the adiabatic olA, When the cycle is 
completed, the area described is the infinitely long 
figure ABfioooA. Since the only heat supplied is r 
during the expansion AB, and no heat whatever is 
rejected, all the heat supplied is converted into work. 
The only way to convert into mechanical work all the heat 
communicated to a substance is to cool the body to abso- 
lute zero temperature. If the cycle be stopped some- 
where before, for instance at ha, some heat must be 
rejected, and the work is reduced to the area ABba. 
According to the definition of the absolute scale of temper- 
ature, equal areas, such as abVa\ a^W^a", etc., corre- 
spond to equal intervals of temperature. Hence, the 
area ABCD is to the total area ABpooaA as the 
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range of temperature dT is to the total range T. There- 
fore 

Audp ^dT 

r T' 

from which Eq. (5) follows directly. 

Prob. 2. — Show that the " physical dimension " on the right- 
hand side of Eq. (7) is " volume per unit weight," and consequently 
checks with the physical dimension of u, Solviion. — The dimen- 
sion of r is — r-TT, that of ^4 is , . Substituting, we get 



weight 



u = 



heat 
weight 



work* 



work 



heat . pre^ure weight X pressure 

, X temp. X . 
work temp. 

_ weight X length _ (length)^ 

' "*' - (^ ' ""*' ■ 

Prob. 3. — For a certain substance the total heat of vaporiza- 
tion may be represented by an empirical formula 

where/ and g are constants; the relation between the pressure and 
the temperature is given by Rankine's formula (20) of the preceding 
chapter. Express w as a function of the absolute temperature T. 
Ans, u = (/log T - gT)/A log-^ (A - B/T - C/T^). 

Prob. 4. — Using the data of Prob. 1 determine the rate at which 
the pressure of steam varies with the temperature at and about 
150° C. Ans, 1306 kg. per sq. m., per degree C. 

Prob. 5. — Pressures p and heat of vaporization r of saturated 
steam at and near 180° C. are given in the following table. 



°c. 


p.kg.persq.cm. 


r in kg-cal. 


175 
180 
185 


9.133 
10.260 
11.493 


482.66 
479.00 
475.34 



From these data calculate the weight of one cu. m. of steam at 
180° C. One cu. m. of water at this temperature weighs about 
887.5 kg. Ans. 5.227 kg. per cu. m. 
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Empirical Formula for r. 

The amount of heat r necessary to vaporize unit 
weight of water without change in temperature may 
be quite closely represented by the following empirical 
formula : 

r = (a - 60 - {cV + d<8) (8) 

Here a, b, c, and d are empirical coefficients, and t is 
the temperature. The first two terms are the prin- 
cipal ones, the last two being comparatively small 
correction terms. It will be thus seen that r decreases 
with increasing temperature. This means that it takes 
less heat to convert one kilogram of water at 200° C. 
into steam than it does at 100° C. The following table 
gives an idea of this decrease of r with temperature: 



t in degrees 
Centigrade. 


r in kg-cal. per 

kg. weight of 

steam. 


100 
130 
160 
200 


536.50 
515.15 
493.56 
464.30 



The values of the four coefficients in Eq. (8) may be 
determined by writmg down this equation for four 
sets of data given in the table above. This gives four 
simultaneous equations of the form 

r2 = a — 6<2 — cfe^ — dU^y 
n = a — htz — ctz^ — dtz^, 
U = a — bU — cU^ — dt^. 



(9) 



To solve these equations in the ordinary way leads to 
wearisome computations of doubtful accuracy. With 
the temperature varying between 100 and 200 degrees, 
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the value of ^ entering into the last terms has an order 
of magnitude of millions; but as the term itself is a 
small one, d must be exceedingly small. Therefore it 
is convenient to represent formula (8) in this form: 

r = a-bt-C {t/l(X)y - D (V100)^ . (10) 

where both C and D are now constants of reasonable 
magnitude, and {t/lOOy varies between 1 and 8. 

The next step is to determme some approximate 
values of a and 6, for instance as in Prob. 6 below. 
Let these approximate values be a' and 6', and let the 
correct values be 



a = a' + a 
b = b' + p 



j, (11) 



where a and /3 are small corrections. The approximate 

formula for r is then 

r' = a'-fe'« (12) 

Subtracting Eq. (12) from Eq. (10) we obtain 

r-r' = a-pt-C (t/lOOy - D (t/lOOy. (13) 

Appljong this equation to the four given sets of values 
gives 

a-pti-C (<i/100)« - D(ti/100y = n -n', \ 
a- pti-C {ti/lOOy - D (<2/100)» = rj - n', 
a-fiU-C (tt/lOO)* - D (<j/100)» = r, - r,', 
a-pU-C (tt/lQOy - D (<4/100)» = u- u'. 



(14) 



In these equations all the fs and r's are given in the 
table, while the values of r' can be calculated from 
Eq. (12). Comparing Eqs. (14) with Eqs. (9) we notice 
that Eqs. (14) need not be solved as accurately as Eqs. 
(9), because all the quantities, viz., a, j8, C, D, are 
correction terms, and a considerable error in their 
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magnitude affects tlie numerical value of r compara- 
tively little. Moreover, all the terms in Eqs. (14) are 
approximately of the same degree of magnitude and 
therefore can be handled more readily and with less 
chance of mistake than the heterogeneous terms in 
Eqs. (9). 

This method for determining the numerical values 
of empirical coefficients has a wide field of appUcation 
in physics and in engineering research, and the student 
will do well to make its principle clear to himself. 

Piob. 6. — Plot the above-given data for r on a sheet of cross- 
section paper and determine the equation of the straight line that 
best fits these points. Ans. a' = 609.1; 6' = 0.724. 

Prob. 7. — Using the equation of the straight line determined in 
Prob. 6 solve Eqs. (14) and determine the coefficients in Eq. (10). 
Ans, a = 606.2; b = 0.688; C = 0.865; D = 0.12. 

Prob. 8. — In order to simplify the solution of Eqs. (14) the 
coefficient Z> is at first assumed to be equal to 0, and approximate 
values of «', /3', C are calculated. Show how to determine the 
accurate values of a, fi, C, and Z), using these approximate values. 

Absolute Error and Relative Error. 

Let it be required to Ifind the maximimi error due 
to using formula (12) instead of (10). In Fig. 16 the 
heavily drawn curve represents Eq. (10), and the slanted 
straight line which crosses it at two points represents 
Eq. (12). It will be seen from the figure that the maxi- 
mum error occurs either at the ends of the range of 
temperatures or at the point where the tangent to curve 
(10) is parallel to straight line (12). The difference of 
ordinates at the ends of the range can be calculated 
directly from Eqs. (10) and (12). The difference at 
point N can also be calculated, provided that the 
abscissa of this point is known. 
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(15) 



The slope of curve (10) at any point is 

^= _!,_ 2(7/ t \ ZD / t \* 

dt loouooj looviooj' ' 

At point N this slope must be equal to the slope —b' 
of the straight line (12). We thus have 

_fe_2C/_fc_\_3^/jLV= -h' 

lOOUOOy 100 UOO/ ' * 

a quadratic equation from which the abscissa to may be 
calculated. Then the diflference r — r' can be calcu- 



(16) 




Temperature t 

Fig. 16. — An accurate and an approximate curve for the heat of 

vaporization. 

lated from Eqs. (10) and (12). The value of this 
difiference giye^ the absolute error mvolved in using 
the simpler expression for r. The relative error is 
equal to (r — r^)/r. 

Prob. 9. — Show that Eq. (16) may be written directly, by mak- 
ing the expression for the absolute error (r — r') a maximum, in 
other words by equating to zero the first derivative of (r — r'). 

Prob. 10. — Show how to find the maximum relative error between 
Eqs. (10) and (12). Hint, (r - r')/r = 1 - (r'/r), so that it is 
sufficient to find a minimum of the expression r'/r. 
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Critical Temperature. 

Expression (10) and the table of values of r show that 
r decreases when the temperature increases. At a 
certain temperature r becomes equal to zero, and the 
distinction between the liquid and its vapor vanishes, 
A very small addition of heat converts the Uquid into 
gas and vice versa. Above this temperature, known 
as the critical temperature, the substance is supposed 
to exist in gaseous state only, and theoretically could not 
be condensed by any pressure whatever. Some experi- 
ments at extra high pressures throw considerable doubt 
upon the last conclusion; nevertheless, the concept of 
critical temperature has proved to be very convenient 
and fruitful in theoretical thermodynamics. 

Prob. II. — Assuming the general form of Eq. (10) to hold true 
up to the critical temperature, show how to determine this temper- 
ature by the method of successive approximations. Perform the 
actual numerical cpmputationa, assuming r per kilogram of water 
vapor to be given by the expression 

r = 606.5 - 0.695 «- 0.2 (^)'- 0.3 (4)'. 
where t is the temperature on the Centigrade scale. Ana. 720° C. 



CHAPTER V. 
ENTROPY OF GASES, LIQUIDS, AND VAPORS. 

The state of a body is determined thermodynami- 
cally by its volume, pressure, temperature, and intrinsic 
energy. The last characteristic is the energy possessed 
by the molecules of the body over and above that pos- 
sessed by the same body in another state called the refer- 
ence state at which this energy is arbitrarily assimied to 
be equal to zero. In addition to these foiu* characteris- 
tics a body possesses a fifth cardinal property, called its 
entropy. This quantity is practically unknown in our 
everyly experience, but tt pUy. L imporUn. port, 
in the analysis of heat engines. 

Perhaps the simplest way to explain the physical 
meaning of entropy is to assiune heat phenomena 
to be due to a hypothetical heat-fluid. Let this fluid 
permeate all the bodies of the universe, as a sponge 
is soaked in water. This hypothetical fluid may be at 
a different pressing, and this pressure we conceive of 
as the temperatiu'e of the body. The quantity of the 
hypothetical heat-fluid in a body is its entropy. The 
quantity of heat energy in a body, according to this way 
of looking upon heat phenomena, is the energy possessed 
by its heat-fluid " compressed " to a certain temper- 
ature. To conmiunicate a quantity of heat to a body 
means to add to it a certain amount of entropy brought 
to a certain temperature. This entropy either in- 
creases the intrinsic energy of the body or is utilized 

78 
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in performing external work. In the latter case the 
common expression that " heat disappears, being con- 
verted into work " becomes: " the entropy loses its 
pressure (temperature) and its potential energy de- 
creases, being converted into external work." See the 
note at the end of the chapter. 

Mathematically, entropy is defined by the expression 

dH = T'd4>y (1) 

where dH is an infinitesimal quantity of heat added to 
the body, d<f> is the corresponding increase in entropy, 
and T is the absolute temperature of the body at which 
heat is added. 

Expression (1) is consistent with the foregoing ex- 
planation of entropy and may be made clearer by 
applying it to other branches of physics; for instance, 
the amoimt of energy dH necessary to add a quantity 
of water d<t> to a tank whose level is T feet above that 
of the supply reservoir is T d<f>; the amoimt of work 
dH necessary to increase the electric charge of a con- 
ductor by an amount d<t>y the conductor being at a 
potential T, is T d<t>. In general, energy may be re- 
solved into a qaardity factor, d</>, and an intensity factor, 
T. Eq. (1) expresses the fact that heat energy may also 
be resolved into two such components. The tempera- 
ture of the body is the intensity factor and the increment 
of entropy, 

d« = f (2) 

is the quantity factor. An infinitesimal increment of 
entropy, d<^, must be used in the foregoing definition, 
and not the total entropy, <^, because the temperature, 
generally speaking, changes during the addition of heat. 
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Similarly, when adding water to a tank, its level is 
raised and it is necessary to write an equation of 
work done during an infinitesimal addition of water, 
during which th^ level may be assimied to remain 
constant. Total energy, H, is foimd by integrating 
Eq. (1) between the required limits, and total entropy, 
(t>, is found by integrating Eq. (2).* 

If entropy be considered as the quantity of a hy- 
pothetical heat-substance in a body, then entropy is 
one of the characteristics of the state of the body, the 
same as its volume, pressure, or temperature, for it 
is impossible to add entropy or heat-fluid to the body 
without changing either its temperature, volume, or 
pressure. Moreover, the difference in the amount of 
entropy between two states is a definite quantity, de- 
pending only upon the states themselves, and not on the 
way or path by which the body was brought from one 
state to the other. This is explained in Fig. 17. Let 
A and B represent two states of a body, referred to 
pressure and volmne as coordinates. It is possible to 
change the conditions from A to B in an innumerable 
number of ways; curves 1 and 2 represent two such 
paths. Let the ordinates of curve aft represent the 
heat added along path 1, so that kPi represents the 
total heat added in passing from A to -B ; let curve aft 
have a similar significance for path 2. Let aa repre- 
sent the entropy in state 1. In order to calculate the 

* The treatment in this chapter is limited to so-called reversible 
changes, that is, to slow, gradual changes of the state of a body; more- 
over the changes are assumed to be such that they can be performed in 
the reverse direction, the result being to reduce everjrthing inside and 
outside the body to the initial conditions of heat, work, temperature, 
volume, pressure, position, etc. A consideration of irreversible changes, 
while of extreme importance in physics and in heat-power engineering, 
is outside the province of this work. 
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entropy kb in the state B, it is necessary, according to 
Eq. (2), to take the sum of infinitesimal increments dH 
divided each by the corresponding temperature (isother- 
mal curves are shown in the sketch by dotted lines). 



1>,H,0. 



N \ \ iv "- N » ^s 




— 1> 



Volume 

Fig. 17. — Illustrating that a change in entropy between two states 
is independent of the external work done. 

In so doing, two different curves 6b are obtained, 
but the final ordinate A;6 is the same in both cases, so 
that the difference M of the entropy in the two states 
is a definite quantity independent of the path by which 
the change is accomplished. 

On the other hand, the quantity of heat added is a 
fimction of the path, because this heat not only in- 
creases the intrinsic energy of the body, but also per- 
forms external work, represented by the area oABk. 
This work is different for the two paths, hence, according 
to the law of the conservation of energy, the heat added 
must also be different. However, subtracting the work 
performed from the total heat added, the difference is 



82 ENGINEERINQ MATHEMATICS, [Chap. V. 

independent of the path, and the result gives two curves 
aCf the increase kC of the intrinsic energy at B over 
that in state A being a definite quantity. The im- 
portance of the concept of entropy lies largely in the 
fact that it permits of dealing with changes in heat 
content without considering the external work done. 

It is shown below how to calculate the increase in 
entropy over that in an assumed reference state, in 
application to perfect gases, to liquids and to saturated 
vapors. Printed tables of entropy of the most impor- 
tant Uquids and their vapors are available, and are used 
in the design and the analysis of performance of steam 
engines, steam turbines, refrigerating machinery, etc. 

The method pursued in computing entropy is that of 
Eq. (2), in other words an expression is found for dH, 
then it is divided by the corresponding T, and the result 
is integrated between the desired Umits. We thus 
have for a change in entropy between two states 

The subscripts 1 and 2 refer to the initial and the final 
states of the body respectively. 

Entropy of a Perfect Gas. 

The expression for dH for a perfect gas is given in 
Eq. (2), Chapter II. Dividing both sides of the equa- 
tion by T we obtain 

, dH dT ^ .pdv 

In order to integrate this equation, the ratio p/T 
must be eliminated from it by means of the character- 
istic Eq. (1) of the gas, in the same chapter. This gives 

d<^ = c„y +Ai2-^, .... (4) 



Chap. V.] ENTROPY OF GASES, LIQUIDS, VAPORS. 83 

or after integration, 

02 — <^i = c„Ln ^ + ARLn—. . . (5) 

il Vi 

If entropy 0i in the initial state, characterized by Ti 
and Vif is given or is assumed equal to 0, entropy ^ in 
any other state may be calculated from Eq. (5). It 
is convenient to select for the reference state one in 
which the entropy is less than in any other state dealt 
with in a certain class of problems. Then in prac- 
tical applications entropy always remains positive. 
The absolute zero temperature could not be selected 
as such an initial state, because the entropy coimted 
from this state would be infinitely large for any other 
state. An intelligent student will find good material 
for thought in this statement. 

If the gas expands without the addition of heat from 
outside, its entropy remains constant. This follows 
either from Eq. (2), or from the physical significance 
of entropy as heat fluid. In this case Eq. (5) gives 

c.Ln'^+ARLn^ =0, 

which agrees with Eq. (13) in Chapter I. This is to 
be expected, since in both cases an adiabatic change is 
considered, or one in which no heat is communicated 
to the gas or abstracted from the gas. 

As has been previously mentioned, the state of a 
perfect gas is determined by two of its characteris- 
tics, for example, the pressure and the volume, be- 
cause the temperature may then be calculated from 
the characteristic equation pv = RT, and the entropy 
from Eq. (5) or (6). Nevertheless, it is sometimes 
convenient to be able to read all these four character- 



84 



ENGINEERING MATHEMATICS. 



[Chap. V. 



istics of a gas simultaneously on the same diagram. 
This may be done, as shown in Figs. 18 and 19, by 
selecting two of the quantities as coordinates, and plot- 
ting the other two quantities in the fonn of a series of 
curves. Thus, in Fig. 18, p and v are selected as co- 
ordinates, while T and are plotted in the form of 




X 

Volume 



Fig. 18. — A diagram in which all the characteristics of a perfect 

gas may be read at once. 

curves. In one set of curves the temperature is con- 
stant for each curve; in the other set the entropy is 
constant for each cm-ve. In Fig. 19, temperature and 
entropy are selected for reference axes, while pressiu'e 
and volume are represented by two sets of curves. 
Having diagrams of this kind it is easy to characterize 
changes of state both quaUtatively and quantitatively. 
For instance, in Fig. 18, the state a is characterized by 
the volume Ox, pressure xa, temperature corresponding 
to isothermal rs and the value of entropy correspond- 
ing to the adiabatic mn. Going to a state 6 the volume 
remains constant, while the pressure, entropy and tem- 
perature increase. In Fig. 19, line km corresponds to 
constant temperatiu'e. Going from k to m the en- 
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tropy and the volume increase, while the pressure 
decreases. 
Isothennals in Fig. 18 are rectangular hyperbolse 



2 
B 

I 








O Entropy 

Fig. 19. — A diagram in which all the characteristics of a perfect 

gas may be read at once. 

whose equation is pv = RT, where T is constant for 
each curve.* 

Prob. I. — Assuming the entropy of air to be zero at the freez- 
ing point of water and at atmospheric pressure, what is its entropy 
at i of the original volume and at a temperature of 150° C? For 
the values of c^, A, and R see Chapter I. 

Ans, <h = —0.0211 imit. 

Prob. 2. — Transform expression (5) so as to have entropy ex- 
pressed through temperatures and pressures, instead of temper- 
atures and volumes. Show that the result is 



<h — <t>i = CpLn^ — ARLn^ , . . 
where Cp = Cv + AR, according to Eq. (10), Chapter I. 



. . (6) 



* The properties of a perfect gas may also be represented by means 
of surfaces in three dimensions, taking, for example, v, 4>, and T for co- 
ordinates. A model of this kind has been constructed by the author 
and is described in Sibley Journal of Engineering, Vol. XXIX, October, 
1914, p. 12. 
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Prob. 3. — Write down expressions for entropy when either 
the temperature, volume or pressure of the gas remains constant 
during the change. Ans. (a) ARLn (vi/vi) + const. 

(6) c„Ln (T2/T1) + const, 
(c) c„Ln {T2/T1) + ARLn {vi/th) 
= CpLn {T2/T1) + const. 

Prob. 4. — Show that for larger values of T' in Fig. 18, isothermals 
are further from the origin. 

Prob. 5. — Show that in Fig. 18 the distances of the vertices 
of the hyperbolae from the origin are as the square roots of the 
temperatures. 

Prob. 6. — Show that when ispthemial curves for equal intervals 
of temperature are plotted to a logarithmic scale the hyperbolae 
become a series of equidistant straight lines inclined at an angle of 
45 degrees to both axes of coordinates. 

The equation of lines of constant entropy or adia- 

batics is 

ptf = const., (7) 

according to Eq. (16) in Chapter I. For the purpose 
here m view, these Unes (Fig. 18) must be drawn so 
as to give an increase in entropy in equal steps. To do 
so the value of the constant in Eq. (7) must be selected 
properly for each curve. This can be done by deduc- 
ing Eq. (7) directly from Eq. (5). Namely, substituting 
in Eq. (5) the value of T = pv/R from the characteris- 
tic equation, we obtain 

^ - 01 = cLn^??^ + AKLa% 

PlVi Vi 

or 

cfe - <Ai = c«Ln^ + (AR + c.) Ln^. 

Pi Vi 

Dividing both sides of the equation by c„ and remember- 
ing that according to Eq. (10) Chap. I, 

Cv + AR=^ Cp, 
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and that the ratio of Cp/Cv is denoted by k, we obtain 

Cv Vi ^i 

Getting rid of the logarithms, the expression becomes 

Pi \Vll 

In this expression pi, Vij and <t>i are the values corre- 
sponding to an arbitrarily chosen reference state. The 
values p2, t^, and <t>2 refer to some other state. If the 
two states are characterized by having the same entropy 
(<fc = 0i), the right-hand side of Eq. (8) becomes equal 
to imity, and the equation itself becomes identical with 
Eq. (7). By giving various values to (fn, which in- 
crease in arithmetical progression, the values of the con- 
stant in Eq. (7) are obtained for the required adiabatic 
curves. In other words, the equation of a family of 
adiabatics may be written in the form 

V2th^ = Vi^i • e^**'*'^^"^, .... (9) 

where the quantities with the subscript 1 are for a 
certain reference state and are constant, while the 
quantities with the subscript 2 refer to any arbitrary 
variable state. For a constant ^ Eq. (9) becomes 
identical with Eq. (7). 

Let, for instance, the reference state for air be that 
at a pressm-e of one metric atmosphere ^and at a tem- 
perature of freezing water. At this state, let the 
entropy 0i per kilogram of air be assumed equal to zero. 
Then for an adiabatic line with entropy ^ = 1 unit, 
Eq. (9) becomes 

p^.405 = 10^000 (0.7733)i-«^ X €i-~/oi», 

specific volumes being in cubic meters per kilogram 
and pressures in kilograms per square meter. Eq. (9) 
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shows that in order to have adiabatic curves corre- 
sponding to equal increments of entropy the constant 
m Eq. (7) must vary in a geometric progression. 

The reference curves shown in Fig. 19 may be an- 
alyzed in a similar manner. At a constant volmne the 
relation between entropy and temperature, according 
to Eq. (5), becomes 

T 
<t>2 — <t>i = c„Ln j=^ + const. . . . (10) 

The shape of these logarithmic or exponential curves 
is shown in Fig. 19. They can be converted into 
parallel straight Unes by using a logarithmic scale for 
T. For t^ = t;i the constant is zero, for t^ > Vi it is 
positive, for t^ < t;i it is negative. 

Prob. 7. — Show that the system of curves corresponding to 
Eq. (9), when plotted to a logarithmic scale, is transformed into a 
series of straight lines inclined to the axis of abscissse at an angle 
whose tangent is — «. Are these straight lines equidistant like the 
isothermals? 

Prob. 8. — On a sheet of logarithmic cross-section paper prepare 
a chart similar to Fig. 18, for one kilogram of air, from the above- 
defined reference state up to a pressure of, say, 50 metric atmos- 
pheres and to a specific volume equal to about 20 times that of 
the reference state. 

Prob. 9. — Replot the chart obtained in the preceding problem 
on a sheet of ordinary coordinate paper. Note. — Whenever a 
curve may be converted into a straight line by plotting it to a log- 
arithmic scale it is usually most convenient to plot it so first, and 
then simply transfer the points to an ordinary scale. It saves 
time and lessens the chance for error. 

Prob. 10. — Using a logarithmic scale for T^ draw lines of constant 
volume for air according to Eq. (10), for temperatures from 0® C. 
to 1000° C. (0° C. corresponds to the freezing point of water), and 
for a range of entropy within one half imit. 

Prob. II. — Show that in the preceding example the distances 
between the lines vary as the logarithms of volumes. On the 
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basis of this relation devise a simple graphical method for draw- 
ing any desired number of these lines, knowing only two points on 
one of the lines or the angle of inclination with respect to the 
T-axis. 

Prob. 12. — Complete the chart obtained in Prob. 10 by draw- 
ing lines of constant pressure according to Eq. (6). Note the 
minus sign before the last term; it means that lines corresponding 
to higher pressures lie above those of lower pressure, while Hues 
corresponding to greater volumes lie below those of smaller volume. 

Prob. 13. — Replot the curves obtained in Probs. 10 and 12 to 
an ordinary scale by transferring ordinates. 

Prob. 14. — Show that for polytropic processes (Chapter II), 
changes in entropy depend only upon the limits of temperature and 
upon the exponent n of the process. 

Entropy of Liquids. 

Roughly speaking, the specific heat of a liquid is 
constant within wide Umits of temperature, so that the 
increase in heat may be assiuned to be proportional to 
the increase m temperature, or 

dH = cdT (11) 

Hence, according to Eq. (2), 



or, integrating, 






•<^-0i=.cLn^^ (12) 



It is customary to count entropy from the freezing 
point of water, so that the above equation becomes 

« = cLn^-ip (13) 

Here <t) denotes the increase in entropy above that 
possessed by one kilogram of liquid at 0** C. 
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The following analogy (Fig. 20) may help to elucidate 
the relation between temperature, entropy, quantity 
of heat, intrinsic energy and specific heat. Entropy 
or heat-fluid may be compared to a ponderable fluid 




Fig. 20. — An analogue for entropy, temperature, quantity of heat, 

and specific heat. 

which is being raised from the level 00' to fill the vessel 
OabO\ The vessel is supposed to have the form of a 
surface of revolution about the axis xy. The temper- 
ature (absolute) is analogous to the difference of level 
between the reference level 00' and the variable level 
mn of the raised hquid. The mechanical work neces- 
sary to raise a quantity d<t> of the Uquid is dH = T 'd<i), 
which expression is identical with Eq. (1). 

The shape of the vessel is determined by condition 
(11), namely, that equal amoimts of work should cor- 
respond to equal increments of level; in other words, 
in proportion as the level of the fluid is raised, the cross 
section of the vessel decreases. Let the cross section 
of the vessel at a height The s; then we have 

d4> ^ 8* dT, 
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or 

dH = cdT = Td<t> = T^sdT, 
thus 

« = f; (14) 

in other words, the cross section of the vessel decreases 
inversely as its height. Eq. (14) gives a simple in- 
terpretation of specific heat c; namely c = sT is the 
volmne of the cylinder pmnq, constructed on the cross 
section of the surface of the vessel. This volmne is 
the same for all cross sections, such as p'm'n'?', because 
of the property of the surface, expressed by Eq. (14). 

Incidentally, this analogue shows why entropy could 
not be coimted from the absolute zero, because at the 
level 00' the cross section of the vessel is infinite, so 
that the total amount of entropy in a body is infinite. 
However, we are concerned only with finite changes in 
entropy. Possibly this analogue could be extended to 
the case when the volmne of the body changes during 
the process, but it would become more compUcated 
and artificial. 

The assumption that the specific heat of a liquid does 
not depend upon its temperature is not accurate enough 
for some purposes, and in calculating entropy it some- 
times becomes necessary to take into account changes 
of specific heat with temperature. Therefore, instead 
of Eq. (12), we have in general 

<h- 4>i = J c—y .... (14a) 

where c is the variable specific heat of the liquid. For 
several Uquids the amount of heat necessary to raise 
their temperature has been represented by empirical 
formulae of the form 

H = ct + aP + pfi, . . . . (15) 
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where a and p axe small correction coefficients which 
account for variations in the value of specific heat. 
If the specific heat were absolutely constant, a and /3 
would be equal to zero, and H would be strictly pro- 
portional to t. Temperature t in Eq. (15) is counted 
from 0° C. and the results are in calories. Differentiat- 
ing Eq. (15) we obtain 

dH = cdt + 2atdt + 3fit^ dt. 

Dividing this expression by T, and remembering that 
dl = dT, we obtain 

d,.f .ef +2a?^«ir+3s21^,,.. (16) 

This expression is easily integrated, and we get an 
expression of the form 

02 - 01 = ALnT + BT + CT^T'^' 

or, denoting by the entropy of the liquid above 0° C, 
= AIjnT/273 + B {T - 273) + C {T^ - 2732). (17) 

For some Hquids, as for example water, the specific 
heat varies rather irregularly with temperature, so 
that the experimental results cannot very well be 
represented by an empirical equation like (15). In 
such cases it is customary to assume the specific heat 
to be constant within certain narrow ranges of temper- 
ature, and use expression (12) for the entropy, adding 
the results. Thus, let the specific heat be constant and 
equal to Cn between the temperatures Ti and T2 ; let it 
be equal to C23 between the temperatures T2 and Tz; etc. 
Then, according to Eq. (12) the total increase in entropy 
between temperatures Ti and T^ is 

0n - 01 = CizLnTs/ri + C2zJjnTz/T2 + • • • 

+ C(n-i)nLnr„/ru.i) (18) 
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Prob. 15. — Assuming the specific heat of water between the 
freezing and the boihng point to be constant and equal to unity, 
plot the temperature-entropy curve (13), and compare the result 
with the data given in published steam tables. Note that the 
curve becomes a straight line when plotted to a logarithmic scale of 
temperature. 

Prob. 16. — Explain the reason why less entropy (or heat fluid) 
need be added to raise the temperature of a kilogram of water from 
99** C. to 100° C. than from 0° C. to r C. 

Prob. 17. — Determine the coefficients Aj B, and C in Eq. (17) 
by differentiating it and comparing with Eq. (16). 

Ans, A = c - 2a X 273 + 3/5 X 273^; 

B = 2a-3/3X 2 X273; C = 3/3/2. 
Prob. 18. — Determine the coefficients A, B, and C in Eq. (17) 
by directly integrating Eq. (16). (See Ans. to Prob. 17.) 

Prob. 19. — From Eq. (15) prove that the average specific heat 
between two given temperatures h and k is 

Cave = C + a («i + ^) + /S \{ti + ^2)^ - ktil 

Prob. 20. — From Eq. (15) prove that the true value of specific 
heat at a temperature / is 

Ctrue = C + 2at + Z^. 

Prob. 21. Peabody gives the following values of the specific heat 
of water: 

Degrees Centigrade. Spedfic Heat. 

Oto 5 1.0072 

5 to 10 1.0044 

10 to 15 1.0016 

15 to 20 1.0000 

20 to 25 0.9984 

25 to 30 0.9948 

30 to 35 0.9954 

35 to 40 0.9982 

40 to 45 1.0000 

From these data calculate the increase in entropy between 0° C. and 
45° C, and compare the result with that obtained directly from 
Eq. (13) in which c is assumed equal to imity. What should be the 
value of c in Eq. (13) in order to give the correct value of entropy 
between the above-given limits of temperature? 

Ana. From data, ^ = 0.1525817; 
From Eq. (13), <A = 0. 1525933 ; 

c = 0.99992. 
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Entropy of Saturated Steam. 

The process of vaporization may be thought of as 
one of widely separating the particles of a fluid at a 
constant temperature. This separation is done against 
the forces of attraction between the particles. The 
energy thus spent in the form of heat is contained in 
the vapor in a latent form, and is called the heat of 
vaporization. According to our conception of en- 
tropy, the process may be mmgined to consist in supply- 
ing excess heat fluid to the Uquid, this heat fluid being 
associated with vaporized particles. The liquid and 
its saturated vapor may be in equilibrium together, in 
which case the vapor is called wet, and the per cent vapor 
in the mixture is called its quahty. Thus, if all the water 
is evaporated, the quality of steam is 100 per cent; if 
20 per cent of the water is left suspended in steam, the 
quality is 80 per cent, etc. 

Let the quality of steam be x and the heat of vapor- 
ization r. Then the heat necessary to convert one 
kilogram of water into steam of quality x at the same 
temperature T \b rx] hence the amoimt of heat fluid 
added, or additional entropy due to vaporization, ac- 
cording to Eq. (2), is 

xr 

The amount of heat fluid or the entropy 0j of a kilo- 
gram of water just before the vaporization begins is 
expressed by Eq. (12), or more generally by Eq. (14a), 
so that the total entropy of the mixture is 

<t>^<t>i + Y (19) 

Values of entropy of the liquid and of heat of vapor- 
ization are foimd in steam tables, so that ^ may be 
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calculated from (19) for any mixture characterized 
by its quality x and temperature T. In some cases 
the pressure may be given instead of the temperature, 
since for a saturated vapor the pressure depends upon 
temperature only, and the relation between the two is 
foimd in steam tables ; see also Fig. 13. 

The question of entropy of steam is not discussed 
any further here, since it affords no appUcation of the 
calculus. Practical problems are mostly solved by 
means of a steam table or a chart plotted from experi- 
mental data. See for example the MoUier diagram 
which may be found in many standard works on the 
subject, and a more complete and more convenient 
chart devised by Prof. F. O. Ellenwood in his " Steam 
Charts '' (Wiley). 

Prob. 22. — Dry steam (quality 100 per cent) at an initial tem- 
perature of 150° C. expands adiabatically imtil its temperature 
drops to 110° C. What is the quality of the steam in the final 
state? From Zeuner's steam tables, calculated from the results 
of Regnault's experiments, we find the following values of r and ^i: 



Temp. "C. 


r 


*i 


150 
110 


500.8 
529.4 


0.442 
0.340 



SdvUon. — Since the expansion is adiabatic, the entropy of the 
mixture remains constant; therefore Eq. (19) gives 

4>i -r jTf — <^i -r J,,, f 

from which quality x" may be calculated. 

The same problem may be solved more quickly with the aid of a 
chart which gives values of entropy for different temperatiu*es and 
different qualities of steam. To the given initial condition of the 
steam corresponds a certain point on the chart. Drawing a line of 
constant entropy through this point to the interception with the 
line of the final temperature, the final quality is read off directly. 

Ana. 03.05 per cent. 
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Note: For a thoroughgoing theoretical and practical discussion of 
entropy the reader is referred to the three books given in the list at the 
end of this volume. The far reaching concept of entropy as the logarithm 
of the probability of a thermod3mamic state is originally due to Ludwig 
Boltzmann. His ideas have been of late further developed in an inspir- 
ing way and appUed to new problems by Max Planck; see for instance 
his "Eight Lectures on Theoretical Physics," published by Colmnbia 
University in 1915. 

For a further development of the idea of entropy as a hsrpothetical 
heat-fluid see the author's article entitled "Entropy and Temperature 
illustrated by Analogs," in the General Electric Review, Vol. 16, 1913, 
p. 7. Also reprinted in the Engineering News, Vol. 69, 1913, p. 406. 



CHAPTER VI. 

HEAT TRANSFER IN BOILERS. 

One of the most important problems in the design 
and operation of a steam boiler is to obtain a maxi- 
mum efficiency of heating surface, in other words, to 
transfer to the water a large part of the heat contained 
in the products of combustion of the fuel. 

For our purposes, a simple steam boiler (Fig. 21) 
may be considered to consist of a cylindrical vessel B, 




Fig. 21. — A simple steam boiler. 

filled with water, with a flue in the center through which 
hot gases pass from the combustion chamber A into 
the chimney C* Steam is taken from the dome F. 

Let ti be the temperature of the gases in the combus- 
tion chamber and <2 their temperature at the chimney 

* The combustion chamber is shown separately in order to simplify 
the mathematical relations. 

97 
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end of the flue. The heat given up by the gases in cool- 
ing from ti to fe is transferred to the water and is used 
in converting it into steam. For a given fuel and 
given conditions of combustion, temperature <i may be 
considered as given; the final temperature fe of the gases 
must be as low as is practicable without unduly increas- 
ing the length of the boiler. The exchange of heat be- 
tween the gases and the water natiu'ally increases with 
the difference of temperature between the two, so that 
more heat is transferred per unit surface of the flue 
nearest to the furnace than per unit surface adjacent 
to the chimney. 

Our problem here is to establish a relationship be- 
tween the dimensions of the flue and the variable 
temperature of the gases along it, so as to be able to 
determine the total amount of heat transferred to the 
boiler. This enables us to judge of the amount of 
steam which a given boiler is capable of producing, 
and also furnishes a criterion for the thermal efficiency 
of the heating surface. 

Let to be the temperature of the water in the boiler; 
this temperature depends upon the desired pressure 
of steam (see Fig. 13, Chapter III), and let t be the 
temperature of the flue gases at a distance x from the 
furnace. The quantity of heat, g, transferred in one 
hom* through one square centimeter of the heating sur- 
face, is a fimction of the difference of temperatures, or 

?=/(<-<o) (1) 

Let W be the weight of gases that pass through the 
flue per hour, and let Cp be their specific heat at con- 
stant pressm^e.* The amount of heat given up by the 

* The pressure of gases in a boiler flue is but a trifle below the 
atmospheric pressure. The specific heats Cp and Cv are treated in Chap- 
ter I. 
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gases during an hour in cooling through an infinitesimal 

change of temperature dt is Wcpdt, and this must be 

equal to the amount of heat absorbed by the water 

through the corresponding element dS of the heating 

surface. Thus 

qdS=-Wcpdt, (2) 

where the minus sign shows that the temperature t 
decreases with the increase of S, the surface being 
counted from the furnace end of the flue. The neces- 
sity for the minus sign becomes evident if one forms 
the derivative dt/dS; namely, from (2), 

dt^ _ _ Q 

dS " Wcp 



(3) 



The derivative is negative by virtue of the physical 
phenomenon itself, the general character of the curve 
connecting t and S being shown in Fig. 22. But the 




Distanoe from Fumaoe End 
Fig. 22. — Decrease in the temperature of gases along the flue. 

expression q/Wcp is positive; therefore its sign in Eq. 
(3) must be changed to minus. 

Older Theory. 

Eq. (2) can be solved only if q is given as a fimction 
of t — to' The earlier writers assumed the transfer 
of heat to be directly proportional to the difference 



■J J 

i J J -" , 
-. J •' -> ^ 
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t — toy and later when the results did not check with the 
experimental data, they assumed q proportional to 
(t — toy. Let us assmne in general 

q = k{t-tor, (4) 

where the coefficient k and the exponent n are to be 
determined from experimental results. Eq. (2) be- 
comes 

k{t-toydS = -Wcpdt, 

or, being integrated, 

k n 



This gives 
kS 



(5) 



. Wcp 
or 

kS 



{n - i)Ut - uy-A 



n-ll{t 



-toy-' (fi-fc)»-T • ^^^ 



This is the law according to which the temperature 
should drop along the flue, provided that Eq. (4) is 
correct. The area S may be replaced by the expression 
TT Dx, where D is the diameter of the flue, so that 
Eq. (6) becomes a Telationship between x and t (Fig. 
22). Eq. (6) holds true for any value of n, except that 
for n = 1 it assumes an indeterminate form 0/0. In 
this case Eq. (5) becomes 

kS_ ^ _ p dt 

Wcp JtJ - <o' 

or 

Tj^ = Ln (<i - <o) - Ln (f - to), 
which may also be written in the form 



V • 
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Eq. (6) may be solved for the variable temperature 
by transferring the last term on the right-hand side 
of the equation to the left-hand side. This gives 

kS (n - 1) 1 _ 1 

or 

r kS in - 1) 1 -]-<"-» ,0^ 

Eq. (7) is solved by getting rid of the logarithm. 
Namely, 

ti — ^0 _ ,. 

where for the sake of brevity, 

kS 



r = 



TFCp 
Solving for t gives 

t-U = {k-' to) t^ (9) 

By definition, the efficiency of a heating surface is 
the ratio of the actual heat given up by the gases to 
the total heat which the gases could give up if they 
could be cooled down to the temperature of the water. 
Or, the efficiency, 

WCp (ti — ti) ti — h /1/^^ 

'^~Wc,{t,-U)''t,--W ' • • ^ ^ 

where fe is the temperature of the gases at the chimney 
end of the flue. Knowing S and W, temperatm-e <2 
may be calculated from Eq. (8) or Eq. (9), and then 
the efficiency determined from Eq. (10). 

Let Q be the total amount of heat necessary to con- 
vert one kilogram of cold water, at the temperature at 
which it is pumped into the boiler, into steam at the 
pressure which is maintained in the boiler. Let I 
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be the weight of steam taken from the boiler per hour. 

Then 

QP = TTcp (<i - fe) (11) 

Knowing P and Q, either TT or fe may be determined 
from this equation. 

Prob. I. — The temperature of the gases in the fire-box of a 
boiler is 1150** C, that at the further end of the flue is 320** C. 
The temperature at the center of the flue is 440** C. What is the 
value of the exponent n, assuming that the heat transfer takes 
place according to Eq. (4)? The temperature of the water in the 
boiler is 190** C. Hint. — The problem is reduced to an equation 
of the form a* + 6* = 2, which may be solved graphically. 

Ans. n = 1.75. 

Prob. 2. — In the preceding example the boiler is furnishing 
steam at the rate of 15 kg. per sq. meter of heating surface per 
hour. What is the value of A; in Eq. (4)? Hint. — The value of. 
S/Wcp in Eq. (6) may be determined from Eq. (11), and the value of 
n is known from the preceding problem. Assume Q = 475 kg-cal. 

Ans. k = 0.234. 

Prob. 3. — Show that Eq. (7) may be obtained from Eq. (6) 
by disclosing the indeterminate form 0/0 according to the rule given 
in the differential calculus. 

Perry's Theory. 

Accurate experiments on the performance of heat- 
ing surfaces are few, so that a complete theory of heat 
transfer in boilers cannot be estabUshed at present. 
The available material indicates, however, that neither 
Eq. (6) nor Eq. (7) represents the phenomena with 
sufficient accuracy, so that the assumption that the 
transfer of heat depends only upon the difference of 
temperatures seems to be inaccurate. The most inter- 
esting fact in this connection is that by " forcing " the 
boiler, tKat is by burning more coal, an ahnost pro- 
portional increase in steaming capacity is achieved. 
Thus, within reasonable limits of forcing, the final 
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temperature of the gases and the efficiency remain 
nearly the same. This observation led some engmeers 
to investigate more closely the physical phenomena of 
heat transfer, particularly with regard to the influence 
of the velocity of the flue gases. If a double volume 
of gases can be made to give up twice the amount of 
heat within the same temperature range, the inference 
may be that a double velocity of particles makes the 
heat transfer twice as easy, with the same difference of 
temperatm-e between the flue gases and the water. 

Professor John Perry's theory given below takes the 
velocity of gases into account; a possible generaUza- 
tion of this theory is then suggested as a more ad- 
vanced exercise in mathematical theory. 

Perry assumes that the heat q transmitted to the 
water m one hour, through one square centuneter of 
flue surface, may be expressed as follows: 

q = k8v(t-to)y (12) 

where ifc is a constant, 8 is the specific weight of the gases, 
V is their velocity, and {t — to) is, as before, the differ- 
ence of temperature between the hot gases and the 
water in the boiler. The factor 6 appears because the 
number of molecules hnpinging upon the wall and 
givmg up their heat is proportional to the density of 
the gas. The factor t; is an empirical one; experiment 
shows that a gas gives up its heat more readily when in 
motion, probably on account of the greater agihty of 
the molecules, and the greater ease with which cooled- 
off portions of the gas are replaced by new. 

The same weight of gas passes through all the cross 
sections of the flue, so that the velocity gradually de- 
creases towards the chimney, because the gas becomes 
cooler, consequently denser, and occupies less volume. 
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The velocity therefore varies inversely as the den- 
sity and the product to in Eq. (12) is constant. Mul- 
tiplying and dividing Eq. (12) by the cross section A 
of the flue gives 

But vA is the volume of the gases, and this multiplied 
by the density b gives the weight TT.* Hence, finally, 

q^k^it-U) (13) 

Eq. (13) may at first sight be deemed to be identical 
with Eq. (4) for the case n = 1, but there is a radical 
difference between the two, since in Eq. (4) the transfer 
of heat is assumed to depend upon the difference of 
temperatures only, while in the latter case it is assumed 
to be proportional to the quantity of the flue gases, in 
other words upon their velocity. 

Substituting the value of q from Eq. (13) into Eq. 

(2) gives 

W 

kj{t-to)dS^ -Wcpdt, 

or, after integration, an expression similar to (6) : 

kS 



Ac, 



'^(!iH) ("> 



The fundamental difference between Eqs. (6) and 
(14) is that in one case the curve of drop of temperature 
(Fig. 22) depends upon the quantity W of the passing 
gases, while in the other case it is independent of W. 
Practically, it means that, if Perry's formula (13) is 

* While W represents the weight of gases per hour, v may be ex- 
pressed in meters per minute, or per second, the numerical factor 60 
or 3600 being included in the coefficient k in formula (13). 
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true, the efficiency i? of a boiler, expressed by Eq. (10), 
is independent of the amount of coal or oil bimied. 
In other words, the steaming capacity of a boiler may 
be increased indefinitely without impairing its econ- 
omy. Some experiments seem to indicate that this 
conclusion holds true within reasonably wide limits 
and promises a great saving in space and in the first 
investment for large power houses. 

Prob. 4. — The boiler mentioned in Probs. 1 and 2 is now 
"forced," that is, more coal is burned, so as to get 50 per cent more 
steam at the same pressure and temperature. The temperature 
of the gases at the furnace end being the same as before, what will 
be the final temperature and the efficiency of the boiler according 
to the assumptions (a) n = 1, and (b) n = 2? 

Ans. (a) U = 1145** C; 17 = 0.52 per cent. 
(6) <2 = 236.5*^ C; ij = 95.1 per cent. 



Suggested Modification of Perry's Theory. 

Perry's formula (12) implies that the exchange of 
heat under equal conditions is proportional to the veloc- 
ity of gases. This is evidently not true at the lower 
limit, when v = 0. As long as there is a difference in 
temperature t — to there will be some exchange of heat 
between the gas and the water, even if the gases are 
stationary in the flue. The existence of the velocity v 
makes the exchange more rapid. Therefore Eq. (12) 
might be modified to 

q = k8(v + Vo)(t-to)y . . . (15) 

Vo being an empirical constant. When v = the ex- 
^ change of heat is reduced to its lowest limit 

5o = kdvo {t - to), 
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instead of to zero. The introduction of an extra co- 
efficient vo has another practical advantage, in that this 
constant may be so selected as to represent the results 
of experiments to a greater degree of accuracy. In 
other words, formula (15) is more flexible than (12). 
Formula (15) may be written in the form 

q=-k(dv + Vo8) (f - <o) = fcQ + Vo8\(t - <o). 
The density 5 is a function of temperature, 

where do is the density at the temperature of freezing 
water, and t is the reciprocal of the coefficient of ex- 
pansion of the gas. The reciprocal is used because its 
physical dimension is that of temperature, and the for- 
mulae which follow are more symmetrical. The higher 
the temperature the more the gas expands (at con- 
stant pressm-e), and the smaller its density becomes. 

Substituting the value of 5 from Eq. (16) into the 
expression for g, we get 

W 



«-*(xn^J<'-«> 



or 

where, for the sake of abbreviation, a notation is in- 
troduced 

VodoA\ 



t 



m 



= (l+^)- (18) 



Eq. (17) is a generaUzed form of Eq. (13) and is re- 
duced to it by putting the velocity Vq = 0, or tm = r. 



J 
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Substituting the value of q from Eq. (17) into the 
fundamental Eq. (2), we obtain 

M = _ r ^^+<)^^ (19) 

Acp Jt, it+t„){t-to)' ' ' ^ ^ 
or, after integration, by resolving into partial fractions, 

kS tm — T J tl + tm I T + to ^ ti — to /^/^N 

T~ ^ 4 — rr ^Q ^ , ^ + : — rrLnr — r- • (20) 

ACp tn, + to t + tm tm+to t - to 

When vo = 0, or tm = r, formula (20) is reduced to 
(14). Eq. (20) can be solved for t by trial only, that 
is, by assuming various values of t until the required S 
is obtained. 

Prob. 5. — For the boiler described in Probs. 1 and 2, plot 
curves of distribution of temperatiu*e, assuming in succession n = 1 
and n = 2 in formula (4) ; also plot a curve on the basis of Perry's 
assumption (12). Only the temperatures at the extreme ends of 
the flue are given. What data are lacking in order to enable a 
similar curve to be plotted according to the assumption (15)? 
Suppose the temperature at the middle of the flue (440° C.) to 
be given; can the curve of temperatures be plotted according to 
Eq. (15)? 

Prob. 6. — Check the integration leading to Eq. (20). 

Prob. 7. — The boiler described in Probs. 1 and 2 was fired so as 
to give 22 kg. of steam per hour per square meter of the heating 
surface. The final temperature of the flue gases was found to be 
430° C, that at the center of the flue to be 520° C. With these 
data and the data given in Probs. 1 and 2, calculate the values 
of Vo and k to be used in Eq. (15). The diameter of the flue is 0.5 m., 
its length is 6 m. Take the following mean values: the coeffi- 
cient of expansion a = 1/r = 0.0037 per degree centigrade; Cp = 
0.240; 5o = 1-25. Hint. — In solving Eq. (20) for vo assume in 
the first approximation vo = 0, where it comes under the logarithm 
sign. 

Am, First approximation, vo = 218; k — 5.25. 
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Further Development of the Theory. 

An objection may be raised against formula (15), in 
that it gives q = oc when v = oo. It is more reason- 
able to assume that when the velocity of gases increases 
without limit, the quantity of heat transferred ap- 
proaches a certain finite limit, so that it is wasteful to 
increase the quantity of flue gases beyond a certain 
limit. A simple way to take this into account is to 
introduce in Eq. (15) the factor 

Vl 
V + Vi' 

where Vi is an empirical coefficient. Thus Eq. (15) 
becomes 

q = kB^^^^it-to). . . . (21) 

{v + Vl) 

For small velocities v is small as compared to Vi, so 
that the newly added factor does not appreciably 
modify formula (15). On the contrary, for large 
velocities the ratio Vi/{v + Vi) becomes less and less, 
thus reducing the value of q. In the limit, when 
t; = 00, 



/ v + Vo \ 



so that 



V 



= 1, 



g„-oo = Vik8 (t - to)' 

This expression is only Vi/vo times as large as when 
V = 0. Using expression (21) in Eq. (2) gives 



kS 



— tl — t (tn — tm) (tm — t) j U + tm 



ACj, U-r (tn-r) (tm+to) t + U 

, (tn + to) (r + to) J ti - to (o<y\ 
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where, for the sake of abbreviation, 

tm = ii + -^^^jT, (18) 

and 

«»=(l+*^)r ..... (23) 

Eq. (22) can be solved for t by trial only; in practice 
various reasonable values of t or of the final temperature 
U would be assumed and the corresponding values of S 
calculated from (22) . By this method it is not necessary 
to solve the equation for t. The efficiency of the heat- 
ing surface is always determined from Eq. (10), and 
the weight W of hot gases from Eq. (11), whatever 
formula is used for the determination of the relation 
between S and fe, viz., Eq. (6), (7), (14), (20), or (22). 

Prob. 8. — As the velocity of the gases increases indefinitely, 
show that the heat transmitted, according to Eq. (21), approaches 
asymptotically a finite value. 

Prob. Q. — Check the integration leading to Eq. (22). What ex- 
perimental data are necessary in order to solve this equation for 
the constants k, vo, and Vi? 

Prof. Breckenridge's Correction for the Effect of Soot 

and Scale. 

In all the formulae above given for q, namely (1), (4), 
(12), (15), and (21), the amoimt of heat transferred 
through the wall of the flue is assumed to be a fimction 
of the diflference between the temperature of the gases 
and that of the water in the boiler. This is true only 
when the metal wall of the flue is clean on both of its 
surfaces, there being no soot coating due to imperfect 
combustion, and no scale formed from impiu'ilies in 
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water. The heat conductivity of a thin metal wall is 
so high that its temperature on the side of the gases 
is practically that of the boiling water on the other side, 
the temperature gradient across the wall being ex- 
ceedingly small. Therefore, with the clean flue, it is 
practically correct to assume that to, the temperature 
of the water, is equal to the temperatiu'e of the wall 
where the flue gases come in contact with it. 

But when the flue has a considerable coating of soot 
and scale, both of which are poor conductors of heat, 
there must be a considerable temperature drop across 
the wall in order to carry the heat through the soot and 
the scale. Let t^ be the temperature of the free surface 
of the soot which comes in contact with the hot gases. 
Then g is a function of (t — t') and not of {t — to)] 
t' must be substituted for U in the formulae for g given 
above. 

The quantity of heat which flows across the wall is 
proportional to the difference of temperatures t' and 
Uy and therefore is expressed by the formula 

? = A;i (^' - <o), (24) 

where the coefficient A^i depends upon the thickness and 
the character of the soot and scale deposits. The value 
of g in Eq. (24) is the same as in the expressions for the 
heat exchange between the hot gases and the surface 
of the soot, because the same quantity of heat is con- 
ducted across the wall as is taken from the gases. Thus, 
V can be eliminated from the equations, provided that 
the coefficient fci is known. Take, for instance. Perry's 
expression (13), which now must be written in the form 

S = fcf«-0 (25) 
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Substituting the value of t' from Eq. (24) gives 

W 



or, solving for q, 



« = ^z(^"f;-4 



5 .= 



W 

TTTir 



(26) 



Comparing Eq. (26) with Eq. (13) it will be seen that the 
only difference is in the second term of the denomina- 
tor. The quantity of heat q is reduced according to the 
value of the coefficient ki. For ideal conditions, i.e., 
clean surface, the conductivity ki is infinite, and Eqs. (13) 
and (26) become identical. For a wall that is absolutely 
non-conducting, A^i = 0, and consequently g = 0. Us- 
ing Eq. (26) in Eq. (2) a new expression for the temper- 




Veloclty of GKifies 



Fig. 23. — Possible relation between velocity of flue gases and 

heat exchange. 

ature drop along the flue may be developed in a manner 
similar to that in which expression (14) has been de- 
duced. Correction (24) may be applied to any of the 
expressions for q given in this chapter. 
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The results of an experimental study of the trans- 
mission of heat through the metal of the heating plates 
of boilers are presented in Technical Paper 114 of the 
Bureau of Mines, as compiled by Henry Kreisinger and 
J. F. Barkley. A large niunber of temperature deter- 
minations were made of both surfaces of boiler tubes 
in a Heine water-tube boiler while the imit was in 
operation. 

Prob. 10. — Apply Breckemidge's correction in obtaining Eqs. 
(14), (20), and (22). What additional experimental data are now 
necessary in order to determine the constants? 

Prob. II. — It is physically possible that the maximmn exchange 
of heat takes place not at v = c3o , as is presupposed in formula (21), 
but at an intermediate velocity v = a (Fig. 23), so that the curve 
for q reaches its maximum at t; = a, instead of at v = oo . Suggest a 
modification of formula (21) such that it would give a curve shown 
in Fig. 23. The three characteristic values of q must be 

qv^o = kdvo(t — to); 
qt^a = khV2 {t — to)] 
9osQO = kbVi (t — to). 

Show that with the suggested law for g, Eq. (2) may be integrated 
in finite form, and that the result becomes identical with (22) when 
a = 00 and Vi = Vi, 
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